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SEARCH  PATH  EVALUATION  INCORPORATING 
OBJECT  PLACEMENT  STRUCTURE 


1.  INTRODUCTION 


This  report  describes  a  robust  numerical  method  for  evaluating  the  search  performance  of 
operations  that  involve  a  searcher  looking  for  objects  that  are  possibly  arranged  in  a  spatial 
structure.  The  problem  is  motivated  by  applications  in  mine  hunting,  where  individual  mines  are 
sought  in  the  search,  yet  the  mines  are  believed  to  often  occur  with  a  predisposition  to  certain 
patterns.  The  approach  taken  in  this  development  builds  on  previous  work  in  developing  a 
“theory  of  distributed  search”  for  distributed  sensing  applications  (references  1  and  2). 
Probabilistic  representations  of  spatial  patterns  for  objects  are  combined  with  probabilistic 
representations  of  a  moving  searcher  to  obtain  a  mathematically  simplified  representation  of 
cumulative  search  effectiveness  over  time.  This  context  provides  a  parametric  model  of  search 
performance  in  terms  of  geometrical  placements  and  sensing  parameters.  When  applied  to  a 
mine-hunting  context,  the  searcher  is  the  mine  hunter  and  the  objects  are  the  mines.  In 
distributed  sensing  contexts,  such  as  antisubmarine  warfare  (ASW),  the  searcher  is  the  target  of 
interest  and  the  objects  are  the  distributed  sensors. 

The  objective  here  is  to  construct  a  general  search  evaluation  capability  that  is  independent 
of  the  search  path  generation  process  and  that  can  be  used  to  evaluate  the  relative  performance  of 
alternate  search  path  planning  approaches.  As  such,  the  mathematical  models  developed  must 
encapsulate  the  underlying  critical  dependencies  encountered  in  searching  for  mines,  and  little 
more.  Spatially  dependent  sensing  performance,  spatially  dependent  object  likelihoods,  and 
overlapping  search  patterns  are  all  examples  of  the  non-trivial  dependencies  that  are  critical  to 
the  mine-hunting  application.  By  establishing  these  components  on  a  mathematical  probability 
and  numerical  analysis  basis  (as  opposed  to  reliance  on  application-specific  detail),  it  is 
anticipated  that  the  dual  problem  of  distributed  sensing  will  be  a  simple  extension  of  the  search 
evaluation  capability.  The  application-specific  components  of  mine  hunting  are  then  left  to  a 
preprocess  that  translates  them  into  the  geometry,  dynamics,  and  mathematical  probability 
descriptions  that  are  used  in  the  evaluation  capability  itself. 


1.1  SEARCH  THEORY  IN  HISTORICAL  APPLICATIONS 

The  search  for  objects  under  motion  given  uncertainty  is  a  classic  problem  in  search  theory 
that  has  been  well  studied  in  the  naval  application  domain  (references  3  -  6).  General  search 
expressions  for  the  probability  of  successfully  completing  searches  in  nominal  environments  are 
well  known  (references  3  and  4),  and  extensions  to  determine  optimal  search  strategies  in  these 
nominal  environments  have  also  been  well  established  (references  5  and  6).  This  rigorous 
theoretical  basis  is  useful  to  understanding  the  basic  structure  of  search  evaluation,  but  becomes 
numerically  untenable  (due  to  many  complex  integrals)  in  all  but  the  most  trivial  of 
environments. 


The  employment  of  search  theories  in  analyses  creates  a  meaningful  mathematical  mapping 
from  the  descriptions  of  sensors,  targets,  and  searchers  to  the  probabilistic  representations  of  the 
expected  outcomes  of  the  search.  By  careful  study  of  the  probabilistic  performance 
representations,  operational  commanders  are  able  to  more  effectively  allocate  resources  to  a 
search,  modify  search  plans,  or  create  expectations  on  the  length  of  time  needed  to  achieve  a 
search  goal.  The  extents  to  which  these  procedures  are  automated  are  the  subject  of  many  on¬ 
going  studies,  and  are  not  the  topic  of  this  report.  However,  for  either  a  human  or  automated 
decision-maker  to  make  confident  search  decisions  requires  a  robust  capability  for  assessing  the 
search  performance.  When  the  scenario  description  is  more  complicated  than  the  nominal 
considerations  found  in  standard  theory,  numerical  methods  must  be  employed  to  effectively 
evaluate  the  search.  The  standard  approach  to  dealing  with  the  added  complexity  is  to  discretize 
the  search  space  to  a  level  of  granularity  where  each  cell  becomes  a  mathematically  nominal 
problem,  and  then  numerically  evaluate  the  combinatorial  problem  over  the  union  of  the  discrete 
cells.  Such  approaches  do  not  scale  well  when  applied  to  the  size  of  complex  military  operations 
and,  thus,  more  computationally  robust  methods  are  required. 


1.2  SEARCH  THEORY  IN  DISTRIBUTED  SENSOR  NETWORK  APPLICATIONS 

The  use  of  distributed  sets  of  independent  detectors  in  network  formulations  provides  an 
extension  to  classical  search  problems.  In  ASW  applications,  the  fields  of  passive  sonobuoy 
sensors  used  to  find  moving  targets  are  an  example  of  such  a  network.  In  this  context,  the 
detection  events  on  the  individual  sensors  are  defined  to  occur  over  intervals  of  time,  and  the 
occurrence  of  multiple  detections  over  the  sensor  network  is  used  to  reduce  the  false  alarm 
probability.  The  set  of  detections  must  occur  in  a  spatio-temporal  pattern  that  is  consistent  with 
expected  target  motion  as  the  target  traverses  the  field  of  sensors.  While  not  a  standard  search 
problem,  it  has  been  shown  (reference  2)  that  an  extension  of  the  notion  of  detection  with  respect 
to  searcher  activity  provides  a  probabilistically  meaningful  description  of  sensor  search 
performance.  These  are  then  combined  using  statistical  sampling  techniques  for  multiple  sensors 
to  accurately  represent  the  total  search  effectiveness.  In  this  search  context  the  target  is  the 
searcher  (the  moving  object)  and  the  sensors  are  the  objects  of  search  (the  fixed  objects).  The 
probability  approach  takes  advantage  of  the  spatial  distribution  of  sensors  as  a  prior  on  the 
sampling  process,  and  the  distribution  of  expected  target  characteristics  provides  a  marginal  on 
the  performance  of  individual  sensors  over  time.  In  an  analogous  manner,  the  mine-hunting 
problem  is  similar,  with  the  spatial  distribution  of  mine  locations  as  a  prior  on  the  sampling 
process  and  the  searcher  characteristics  given  deterministically  (equivalent  to  a  delta  function 
distribution  for  the  marginal). 

In  both  the  ASW  and  mine-hunting  paradigms,  object  detection  is  said  to  occur  with 
probability  Pd  if  it  occurs  anywhere  within  a  predefined  temporal  interval  or  spatial  region 
respectively.  Let  Xs  =  {xj  be  the  set  of  positions  of  Ns  objects  (sensors)  with  detection  radius 

Rd.  LetxT  be  the  initial  position  of  a  moving  searcher  (target)  randomly  placed  within  search 
space  ScR'  and  traveling  at  known  constant  velocity  vt.  Then,  the  searcher  position  at 
arbitrary  time  t  e[t0,t0  +  AT]  for  initial  time  to  and  fixed  interval  AT  can  be  written  as 
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xT(t)  =  xT(t0)  +  (t-t0)-vT 


(1) 


Let  Xa  be  the  subset  of  Xs  that  falls  within  radius  Rd  of  the  searcher  (target)  trajectory  established 
over  the  time  interval  AT.  Then, 


xd  =  {  XJ  e  xs|  I XJ  - xt(0  I  ^  Rd’ t  e  (to.' to  +  At)} 


(2) 


Observe  that 

xd  =  { XJ  e  xs|  II XJ  - (  XT (f 0 )  +  (t  - 10) •  vT)||  <  Rd, t  €  (t0, t0  +  At)}, 

=  {xj€Xs|xj€Qt}, 

=  {XJ  eXs||  XT(to)-(XJ-(t-t0)-VT)|<Rd,tG(t0,t0+At)}, 

=  {xJeXs|xr(to)eQJ}. 

where  Qt  and  { Qj  }/j  denote  subsets  of  S  centered  about  the  respective  searcher  (target)  and 

object  (sensor)  relative  motion  trajectories  and  extending  in  any  direction  to  a  distance  of  Rd- 
This  is  shown  graphically  in  figure  1  for  the  case  of  upward  searcher  (target)  motion.  Figure  la 
depicts  the  ensemble  of  these  pill-shaped  regions  Qj  originating  at  the  respective  object  (sensor) 
positions  and  extending  downward.  Figure  lb  depicts  a  searcher  (target)  detection  region  Qt  for 
a  single  searcher  originating  at  the  initial  position  and  extending  a  distance  corresponding  to 
searcher  motion  over  the  AT  time  interval. 


100  Sensors.  Spd  Ratio  0  00,  DT  =  2  00  min 


100  Sensors.  Spd  Ratio  0  00,  DT  =  2  00n 


a.  Target  Placement  Within  a  Sensor  Pill  Ensemble  b.  Object  Placement  Within  a  Target  Pill 

Figure  1.  Detection  Event  Equivalencies 
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In  the  former  case,  search  performance  is  governed  by  computing  the  likelihood  that  a 
searcher  (target)  placement  will  fall  within  the  set  of  sensor  regions.  In  the  latter  case,  the 
evaluation  assesses  the  likelihood  of  object  (sensor)  placement  within  the  searcher  region.  The 
only  difference  between  the  two  cases  is  a  frame  of  reference  for  the  observer;  the  search 
performance  evaluation  in  both  cases  is  identical.  The  latter  case  is  directly  analogous  to  the 
placement  of  mines  within  a  searcher  region,  where  for  mine-hunting  purposes  the  searcher 
region  is  defined  by  the  detection  region  associated  with  the  mine-hunting  platform  search 
trajectory.  The  natural  performance  metric  for  either  representation  is  the  probability  of 
successful  search — that  is,  the  cumulative  probability  of  finding  the  intended  object  within  the 
region  covered  by  the  searcher.  The  underlying  probability  assumption  is  that  object  placement 
is  identically  distributed  (although  not  necessarily  uniform)  for  each  of  the  set  of  objects. 


1.3  SEARCH  THEORY  IN  MINE-HUNTING  APPLICATIONS 

The  mine-hunting  problem  domain  is  characterized  by  several  aspects  that  make  it  distinct 
from  both  the  historical  search  theory  and  the  extensions  to  distributed  search.  First,  when 
considering  the  placement  of  stationary  mines,  the  placement  events  are  persistent.  That  is,  once 
a  mine  is  placed,  it  remains  in  its  initial  location  until  some  form  of  intervention  occurs  to 
remove  or  disable  it.  Even  if  mines  are  properly  identified,  their  presence  remains  in  the  sensor 
observation  space  unless  they  are  removed.  Secondly,  sensor  observations  often  include  the 
detection  of  non-mine  clutter  objects  (false  alarms  due  to  objects  that  are  not  the  object  of  the 
search).  When  such  objects  are  detected,  intervention  may  be  required  to  properly  classify  them, 
and  the  spatial  density  of  this  clutter  over  a  search  region  can  exhibit  spatial  variability. 

The  persistent  placement  event  is  fundamental  to  the  search  evaluation  method  described  in 
this  report.  It  is  consistent  with  the  image  feature  processing  conducted  using  side-scan  sonars. 
For  any  given  sweep  along  a  search  path,  a  composite  image  comprising  individual  scans  is 
produced — see  the  example  shown  in  figure  2,  which  shows  a  typical  side-scan  sonar  mapping  of 
a  swath  of  the  ocean  bottom  and  a  search  path  conducted  to  find  mine-like  objects.  The 
detection  assessments  that  are  performed  are  retained  until  subsequent  sweeps  produce  additional 
features  to  process. 


Figure  2.  Example  of  Search  Path  Employing  Side-Scan  Sonar 
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The  use  of  search  theoretic  methods  to  assess  mine-hunting  performance  for  finding  random 
mines  distributed  uniformly  in  a  region  has  been  well  established  (reference  7).  However,  for 
new  applications,  additional  factors  must  be  considered  in  applying  search  theory  to  mine 
hunting.  The  actual  number  of  mines  in  the  search  region  may  be  unknown.  There  may, 
however,  be  limiting  knowledge  on  the  approximate  number  of  mines  in  the  region  based  on 
intelligence  data,  and  such  data  must  be  able  to  be  used  in  a  probabilistic  manner.  The 
deployment  of  mines  may  not  be  identically  distributed  or  even  random.  Mine  field  structure,  if 
it  exists,  or  any  dependencies  on  deployment  such  as  bathymetric  constraints  drive  the 
interpretation  of  the  expended  search  effort.  Prior  knowledge  of  mine  field  deployment  patterns 
is  especially  useful  when  inferring  mine  likelihood  in  the  region  surrounding  identified  mines. 
Such  patterns  include  barrier-like  mine  lines  or  sets  of  mine  lines. 

The  effects  of  background  clutter  that  can  encumber  the  identification  of  mines  must  be 
taken  into  account.  Figure  3a  illustrates  a  spatially  variable  clutter  density  map  (contours)  and  a 
numeric  instantiation  of  mines  (orange  dots)  and  clutter  (red  and  blue  dots).  Additionally, 
independence  assumptions  on  detection  observations  using  multiple  sensors  or  multiple  passes 
with  the  same  search  asset  must  be  considered.  The  consequence  of  search  path  overlap  must  be 
integrated  into  the  search  evaluator,  as  complex  maneuvering  characteristics  may  exist  in 
systems  relying  on  fusion  from  multiple  contacts  to  reduce  the  impact  of  clutter.  Such  a  search 
path  for  high  density  location  interrogation  is  illustrated  in  figure  3b. 


I  0 


a.  Clutter  and  Mine  Field  Spatial  Characterization  b.  Kinematic  Modeling  of  the  Search  Platform 


Figure  3.  Experimental  Equivalencies 


The  objective  in  construction  of  this  search  path  evaluation  capability  is  to  apply  as  much  of 
the  prior  knowledge  on  mine  field  structure  as  is  practical  and  relevant,  while  retaining  flexibility 
in  representing  vehicle  trajectories  and  sensor  performance.  The  next  section  of  this  report 
describes  the  development  of  the  probability  event  spaces  over  which  search  evaluation 
probabilities  are  calculated.  This  is  followed  by  a  description  of  the  resulting  search  evaluation 
realization  in  a  grid-based  design  for  numerical  calculation.  Results  are  presented  to  articulate 
the  performance  of  the  evaluation  capability  in  terms  of  its  accuracy  in  probabilistic  modeling 
and  its  numerical  efficacy. 
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2.  PROBABILITY  EVENT  SPACES  FOR  SEARCH  PROBLEMS 


This  section  provides  the  theoretical  background  on  which  the  development  of  the  search 
evaluation  methodology  is  based.  A  review  of  necessary  background  on  probability  spaces  and 
spatial  Poisson  processes  is  provided  in  appendix  A.  To  facilitate  the  exposition,  an  initial 
simplified  event  context  is  introduced  and  is  progressively  extended  to  obtain  the  event  space 
realizations  necessary  to  provide  an  evaluation  capability  of  sufficient  utility.  Necessary  detail 
on  the  construction  of  all  the  event  spaces  is  provided  to  illustrate  how  it  is  utilized  within  this 
tool  for  evaluating  mine-hunting  search  path  trajectories. 

2.1  SEARCH  FOR  A  SET  OF  KNOWN  OBJECTS 

A  simple  context  is  considered  first  to  illustrate  the  event  space  construction.  Let  ScR' 
represent  a  geographic  region  in  which  to  conduct  a  search  for  a  set  of  certain  objects.  Assume 
that  the  set  of  objects  is  finite  and  that  each  object  is  known  to  exist  within  this  region.  Let 

M  =  { m,  denote  the  fixed  set  of  objects  placed  within  S,  and  let  Xj  =  (xj,yjy  e  S  be  the 

location  of  each  mj.  Furthermore,  assume  that  these  locations  are  unknown  (hence,  the  need  to 
search)  and  let  each  Xj  be  an  independent  random  variable  placed  in  S  with  probability  density 

function  fj(Xj)  for  Xj  eS,  and  zero  otherwise.  Then,  ^x,,...,xNm^  represents  an  NM-dimensional 
vector  on  SN"  with  probability  density  function 


(4) 


The  independence  assumption  between  the  placement  random  variables  implies  that  each  fj(Xj) 
represents  the  marginal  distribution  of  Xj  on  S.  From  this  construct,  the  event  space  associated 
with  the  placement  of  a  single  object  in  S  is  considered. 

Let  A  c  S  be  a  sensor  scan  region  within  S  in  which  one  attempts  to  detect  the  objects  M. 
The  probability  that  an  object  mj  lies  within  A  is  given  by 


(5) 


Pr(Xj  e  A)  =  |AfJ(xj)dxj  =  PA  . 


Since  the  set  {Xj  e  A,  xf  g  A}  forms  a  partition  of  S  into  the  scan  region  A  and  its  complement 
SYA,  the  probability  that  the  object  is  not  in  A  is  given  by 


Pr(xJ«A)  =  Js/AfJ(xJ)dxj=l-PA. 


(6) 


This  search  experiment  in  A  has  two  possible  outcomes  over  a  decision  space  2)  =  {d,  d} 
corresponding  to  detection  of  the  particular  object  or  non-detection.  Continuing  with  the 
assumption  that  object  mj  lies  within  A,  the  detection  event  =  d  is  defined  to  occur  with 
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probability  Pd  =  Pr(^  =  d  |  X|  e  A)  and  the  non-detection  event  d-  =  d  to  occur  with  probability 
1  -  Pd  -  Pr(d,  =  d  |  e  A).  For  this  context,  object  detection  is  restricted  to  the  scan  region 
and  independence  between  placement  and  detection  outcomes  is  ascertained.  The  joint 
placement-detection  outcomes  (xj,djj  are  elements  of  the  product  space  S  x  5) ,  which  is 

partitioned  as  {Xj  e  A,  \t  g  A}x  {d,d}.  The  event  space  5;  =  a({Xj  e  A,  x^  g  A}x  {d,d})  becomes 
the  a-algebra  generated  by  this  partitioning.  Using  Bayes  theorem  (reference  8),  the  probability 


of  detecting  object  mj  in  A  becomes 

Pr(Xj  €  A,<£j  =  d)  =  Pr(rfj  =  d  |  xs  e  A)Pr(Xj  €  A)  =  PDPA .  (7) 

The  remaining  partition  probabilities  are  calculated  in  similar  fashion: 

Pr(Xj  g  A,dt  =  d)  =  Pr(^  =  d  |  Xj  g  A)Pr(Xj  g  A)  =  0-(l  -  PA)=  0,  (8) 

Pr(x,  g  A ,dy  =  d)  =  Pr (df  =  d  |  \s  g  A)Pr(Xj  g  A)  =  1  (l  -  PA)  =  1  -  PA  ,  (9) 

Pr(Xj  e  A, ^  =d)  =  Pr =  d  |  xs  e  A)Pr(Xj  €  A)  =  (l  -  PD)PA  .  (10) 


Observe  that  non-detection  of  the  object  outside  the  scan  region  is,  by  definition,  a  certain  event 
(as  per  equation  (9))  and  that  the  probabilities  over  all  outcome  partitions  sum  to  unity.  The 
latter  occurs  because  a  partitioning  of  S  x  5)  produces  subsets  that  are  both  mutually  exclusive 
and  collectively  exhaustive. 

The  event  space  is  extended  by  partitioning  S  into  multiple  scan  regions.  Following  the 
analysis  conducted  in  reference  7,  the  search  space  is  divided  into  a  set  of  identical  disjoint 
search  regions  where  each  region  constitutes  an  independent  scan  for  objects.  An  example 
partitioning  is  illustrated  in  figure  4  as  a  set  of  rectangular  search  regions.  This  type  of 
partitioning  might  be  provided  for  covering  a  rectangular  search  space  with  a  side-scan  sonar. 
The  actual  shape  of  these  disjoint  regions  is  arbitrary  as  long  as  their  union  realizes  the  space 
covered  by  a  non-overlapping  search  path  over  which  the  search  effectiveness  is  to  be  evaluated. 


Search 

Space 

S 


Figure  4.  Search  Space  Partitioning 
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Let  { A(  be  an  exhaustive  partition  of  the  search  space  S  into  a  set  of  disjoint  scan 

regions  of  equal  area  a*.  Throughout  the  sequel,  the  area  of  region  Aj  is  represented  as  a* 
according  to 

ai=|A,|=  £  ldx.  (11) 

The  probability  that  object  mj  lies  in  region  Aj,  for  i  =  1,  Ns,  is 

Pr(xJeA,)=  jA  fJ(xJ)dxj  =  PA  ,  (12) 

where  it  is  recalled  that  fj(Xj)  is  the  location  probability  density  function  for  object  mj.  Since  the 
partitioning  is  exhaustive  over  S,  the  placement  of  object  mj  in  S  is  certain  and 

Pr(xIeS)  =  SPr(xjeA,)  =  £p,  =1.  (13) 

i=i  i=i  ' 

For  the  moment,  assume  that  placement  and  detection  outcomes  are  independent  for  each  scan 
region.  Then,  as  above,  the  probability  of  detecting  mj  in  Aj  becomes 

Pr(Xj  e  Aj,^  =d)  =  Pr^  =d  |  \)  e  A^P^Xj  e  A;)  =  PDPA  .  (14) 

For  notational  simplicity,  let  Ej  denote  the  joint  event  Ej  =  {^Xj,^  |  Xj  £  A;,^  =  d}  representing 
the  detection  of  an  object  that  exists  in  region  Aj.  Then,  Pr(Ej)  =  PDPA  and  the  probability  of 
not  detecting  object  m,  in  Aj  becomes  Pr( E ; )  =  1  —  PDPA  ,  where  E;  =  x^lx^A.u^d} 

denotes  the  set  complement  with  respect  to  the  event  partitions.  Note  that  Pr(Ej  n  Ek )  =  0,  i  *  k 
as  random  variables  are  single  valued  in  any  specific  realization  and  the  scans  are  mutually 
exclusive.  However,  also  for  i  ^  k. 


Pr(E|  nEJ  =  P^Ej  uEk)  =  l  —  Pr(Ej  uEJ  =  1-PD(PA  +PA  )• 


(15) 


For  search  path  evaluation,  let  Sp  denote  a  search  plan  comprising  a  subset  of  { A,  .  That  is, 

Sp  =  { A,  }jeIp ,  with  Ip  representing  the  indices  corresponding  to  those  scans  conducted  in  a  partial 
execution  of  the  search  plan. 

Now,  assume  that  any  scanned  region  is  searched  only  once.  Then,  the  probability  of 
detecting  mj  in  one  of  the  set  of  disjoint  scans  comprising  Sp  becomes 


Pr 


(JE,  =  Pr  Xj  e  U Aj  nrf|  =d  =PD  IPA.  . 


vie,p  / 


iel 


ie  1P 


(16) 
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Observe  that  when  Sp  includes  all  the  possible  scans  of  the  search  space  S,  then 


Pr(xJeSnrf|=d)=PDl  =  PD. 


(17) 


If  the  search  is  definitive  in  that  the  entire  space  S  is  searched  and  Pd  =  1,  then  the  object  m,  must 
be  found.  This  is  shown  below  in  the  sequential  update  of  the  probability  of  finding  object  m,  in 
scan  region  k+1,  given  that  it  was  not  found  in  the  previous  k  scans: 


riE, 


i=l 


(~ 

Pr  Ek+1|nE, 

V  i=l  ) 

p{n'E, 


•-PdXPa, 

i=l 


k  _  \ 

Pr|  f|Ej 

0=1  ) 


1-PdZPa, 


PDPA,i 

1- PdXPa, 

i=l 


(18) 


For  this  definitive  search  case,  the  probability  of  detecting  object  m,  in  the  last  scan  region  ANs , 
conditioned  on  not  finding  it  in  any  of  the  other  regions  becomes  a  certain  event,  is 


Pr 


Ns-I _  ^ 

Em,  I  nf 


p  p 

rDrA 


NS 


1PA 


ans 


Nc-1 


i-pdEpa,  1'H,"P*-) 

i=l 


(19) 


which  is  as  expected. 

Next,  to  construct  a  more-general  event  space  for  search  evaluation,  some  of  the 
constraining  assumptions  are  relaxed.  In  particular,  the  requirements  on  the  independence 
between  placement  and  detection,  the  requirements  on  spatial  invariance  of  the  detection 
probability  Pd,  and  the  requirements  on  single-scan  detection  assumptions  are  removed.  To  do 
this,  a  search  space  partition  is  constructed  that  is  a  refinement  of  the  partitioning  according  to 
scan  regions.  Specifically,  the  search  space  is  partitioned  into  a  uniform  grid  over  the  placement 
position  coordinates,  and  probability  cells  are  defined  about  the  grid  points.  Then,  as  above,  the 
probability  that  object  mj  lies  within  grid  cell  (k,l)  is  given  by 


Po„  =Pr(xJeGu)=  J0  fj(xj)dXj  • 


(20) 


Since  the  probability  density  function  f(x,,...,xNy)  is  defined  to  have  support  on  S,  the 
detection  probability  conditioned  on  placement  within  grid  cell  (k,l)  becomes 


10 


(21) 


=  Jcu  Pr(d,  =  d  I  xj)fj(xj  I  xj  e  Gk,l)dXj* 


The  term  Pr(^  =  d  |  x  j )  represents  the  traditional  sensor  performance  prediction  of  the 

probability  of  detection.  Observe  that,  if  sufficient  resolution  in  the  grid  spacing  is  attained  such 
that  fj(x j  |  Xj  e  Gk  l)  is  nearly  constant  over  the  grid  cell,  then 


(22) 


where  |Gk,i|  is  the  area  of  grid  cell  Gy,  as  given  by  equation  (1 1).  Thus,  with  a  fine  grid 
resolution  PDk  reduces  to  a  spatial  average  of  the  detection  probability  over  the  grid  cell. 

Conversely,  if  the  grid  resolution  is  coarse  relative  to  fj(xj),  then  the  dependency  between  the 
detection  and  placement  probabilities  is  imbedded  within  the  calculated  values  in  the  form  of  the 
weighted  average  of  equation  (2 1 ). 

For  PDk  <  1 ,  there  is  a  probability  of  1  -  PD ^  that  object  mj  will  not  be  detected  on  the  first 

scan  opportunity.  It  may,  however,  be  detected  on  subsequent  passes  if  the  search  path  crosses 
(or  overlaps)  its  own  trajectory.  Let  ip  =  {ip },ip  e  {l,...,Ns}  denote  a  scan  sequence  that  may 

include  repeated  scans.  By  modeling  each  scan’s  detection  observation  as  an  independent 
Bernoulli  trial,  the  waiting  time  W  (in  terms  of  numbers  of  scans  before  detection  occurs) 
follows  a  geometric  distribution  (reference  8).  That  is,  the  probability  that  the  first  detection  of 
object  mj  will  occur  in  cell  (k,l)  on  the  nth  scan  becomes 


Pr(W  =  n)  =  PDk,(l  — PDk ,)"'. 


(23) 


The  succession  of  first  occurrence  probabilities  develops  sequentially  as  repeated  scans  are 
included  in  the  search  plan.  Let  d-(n)  =  d  denote  the  occurrence  of  the  first  detection  of  object 

mj  within  scan  n  and  ^(n)  =  d  denote  the  occurrence  of  no  detections  of  object  mj  through  the 
first  n  scans.  With  initial  values 


P^  ,(l)  =  Pr(rfj(l)  =  d  |  k,  1)  =  P, 


(24) 


and 


P^  (1)  =  Pr(^(l)  =  d  |k,l)  =  l-P, 


(25) 


the  scan  progression  updates  for  n  =  2  and  beyond  are  given  by 


P4i(n)=Pr(4(n)  =  d|k,l)  =  PDki.P^(n-l), 


(26) 


and 


(n)  =  FM4(n)  =  d  |  k,l)  =  (1  -  PD„ )  ■  p4  (n  - 1) . 


(27) 


Observe  that  PDk  =  1  for  definitive  search  detection  probabilities.  In  that  case,  all  subsequent 
first  detection  probabilities  in  the  sequence  go  to  zero,  as  expected. 

Without  loss  of  generality,  the  cell  detection  probabilities  can  be  developed  based  solely  on 
a  geographic  basis  to  reflect  dependency  on  the  acoustic  environment.  They  can  also  be 
developed  as  functions  of  a  predefined  vehicular  search  path  to  incorporate  dependencies  on 
proximity  to  the  sensor.  In  more  complex  cases,  any  combination  of  geometric  reference  and 
search  path  for  one  or  more  sensors  may  be  applied.  For  situations  with  specified  vehicular 
search  paths,  the  detection  probability  may  need  to  be  recalculated  for  each  subsequent  detection 
experiment  as  proximity  to  the  sensor  changes — this  is  reflected  in  the  form  of  the  sequential 
update.  For  situations  with  a  strict  geographic  reference,  the  detection  probabilities  need  to  be 
calculated  only  once. 

Given  the  cell  detection  probabilities,  the  event  probabilities  defined  for  the  scan 
partitioning  are  evaluated  over  the  grid  partition.  Let  f  designate  the  indices  of  the  grid  cell 
elements  contained  within  the  scan  corresponding  to  event  Ej.  This  event  may  comprise  any 
component  of  the  search  plan  under  evaluation  and  may  possibly  overlap  other  scans  executed 
within  the  search  plan.  The  probability  of  detecting  object  m,  within  this  scan  is  calculated  by 
aggregating  the  probabilities  over  the  grid  partition  refinement,  as  in 


Pr(E,)=  Z  P,  (n(k,l))-PG 


(k.l)el, 


(28) 


This  calculation  accommodates  scan  overlap  at  the  grid  cell  level.  That  is,  first  detection  occurs 
for  a  specific  grid  cell  if  it  occurs  anywhere  within  that  grid  cell.  The  scan  sequence  probability 
Pv  is  retained  individually  for  each  grid  cell  element  (k,l). 

“k  .1 


The  event  space  definition  for  the  case  of  a  known  number  of  objects  is  completed  by 
considering  the  detection  over  the  finite  set  of  Nm  objects.  An  inference  is  drawn  over  the  set  of 
objects  M  =  { m,  by  constructing  event  spaces  pertinent  to  the  probability  question  under 

consideration.  Consider  the  probability  of  detecting  all  NM  objects  in  S,  which  is  known  as  a 
clearance  probability.  Let 

Ea =n{<xj’<fj)ixieAi’rfj =d!  (29) 

j=i 
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represent  the  joint  event  that  all  objects  lie  within  the  scan  region  Aj  and  are  detected.  This  joint 
event  forms  a  partition  over  SNm  x  5)Nm  whereby  each  S  x  5)  subspace  is  partitioned  identically 
according  to  the  search  plan.  Hence,  unlike  the  previous  calculations  where  only  the  subspace 
under  consideration  was  partitioned,  the  entire  outcome  space  must  now  be  considered  to 
calculate  the  probability  that  event  Eci  will  occur;  specifically, 


Pr(ECj)=nPr(xjeA1,tfJ=d)=n 

j=i  j=i 


I  Pi  (n(k,l))  P^ki 

(k.l)elj 


(30) 


In  this  expression  and  throughout  the  sequel,  the  dependency  of  the  grid  cell  probability 
calculations  on  the  respective  fj(Xj)  is  implied.  To  calculate  the  clearance  probability,  the  search 
region  probability  is  calculated  in  each  subspace  and  combined  in  the  product  as 


Pr(Ec)  =  nPr(xjeUA,^=d)  =  f[ 


j=l 


ielp 


j=l 


X  2X(n(k,l)X, 


ie  lp  ( k,l)el, 


(31) 


This  is  a  consequence  of  the  assumption  of  independence  between  object  placements.  When  the 
respective  fj(xj)  are  identically  distributed. 


Pr(Ec) 


I  X 

ielp  (k,l)el, 


P*,(n(k.l» 


(32) 


Next,  consider  the  event  space  corresponding  to  the  detection  of  one  or  more  objects  in  M. 
Here,  the  same  partitions  formed  over  SNm  x  5)Nm  as  used  for  the  clearance  probability  are  again 
used.  Let 


'Di 


=  U{(xJ,rfJ)|xJ€Ai,^=d} 

j=l 


(33) 


represent  the  joint  event  that  at  least  one  of  the  Nm  objects  lies  within  scan  region  A,  and  is 
detected.  Since  multiple  object  placements  are  defined  to  be  mutually  independent,  set 
complements  are  applied  over  the  disjoint  partitions  of  the  outcome  space  to  get  the  joint 
probability.  In  particular, 


Pr(£Dl)  =  1  -  Pr (Ea)  =  1  -  ri  Prfxj  eA,v  d ,  *  d\ 

7=1 

"»(  \  N"( 

=  l-n(l-Pr(xje  4,4  =«/))=!- n  1- 


z  Pj^  {n(k,l))  P£ 

(k,l)el, 


(34) 


To  obtain  the  probability  of  detecting  at  least  one  of  the  Nm  objects  in  a  partial  search  plan,  scan 
probabilities  are  again  aggregated  and  combined.  As  above,  when  the  fj(Xj)  are  identically 
distributed, 
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(35) 


N«  /  \  f 

Pr(ED)  =  1  -  n Pr  Xjg  UAjUt^^d 

j=l  ^  ielp 


X  X  PA.(n(k,l»-Pc; 


ielp  (k ,l)el, 


\Nm 

J 


as  the  probability  of  detecting  at  least  one  of  the  Nm  objects. 

This  section  concludes  by  focusing  on  the  case  of  identical  placement  distributions,  as  it 
illustrates  an  interesting  property  of  this  event  space.  In  this  case,  the  detection  of  each  of  the 
respective  objects  within  the  search  path  is  an  independent  Bernoulli  trial.  Let  p  be  the 
probability  of  finding  an  object  in  the  search  path: 


P  =  X  X  Pdll(n(k,l)) 

ielp  (k, l)el, 


(36) 


Then,  Pr(Ec)  =  pNM.  Similarly,  with  q  =  1  -  p,  then  Pr(ED)  =  l-qNM  or,  equivalently, 

Pr(ED)  =  qNy .  More  generally,  let  N<j  be  the  number  of  objects  detected  in  the  search  path  and 
ENj=m  be  the  event  in  which  Nj  =  m  occurs.  Assuming  that  the  order  of  object  detection  is 
unimportant,  the  probability  of  event  EN  =m  occurring  is  given  by 


Pr(ENj=m)  =  Pr(Nd  =m)  = 


N 


M 

m 


p'VM-m 


(37) 


Thus,  the  assumption  of  identical  distributions  in  the  placement  of  objects  leads  to  a  binomial 
distribution  on  the  number  of  objects  detected  in  the  search  path.  While  this  event  labeling  has 
changed  to  a  count  indicator,  the  event  space  has  not  changed.  Hence,  Pr(Ec)  =  Pr(Nd  =  NM ) 

and  Pr(  ED)  =  1  -  Pr(Nd  =0).  Note  that  for  sufficiently  large  NM,  as  Nm  — 1 ►  oo  and  p  «  1  such 
that  Xp  =  Nm  •  p  remains  finite,  then 


Pr(EN  )  =  Pr(Nd  =  m) 


ur 


m! 


(38) 


by  the  Poisson  approximation  to  Bernoulli  trials  (reference  9). 

It  is  important  to  note  that  with  identical  placement  distributions,  the  calculation  of  search 
effectiveness  probability  is  effectively  reduced  from  an  operation  on  the  SNm  x  2)Nm  outcome 
space  to  an  apparent  operation  on  S  x  $)  through  a  functional  composition.  That  is,  the  repeated 
operation  on  S  x  ID  is  embedded  in  the  functional  forms  of  equations  (32)  and  (35).  Hence,  the 
mappings  are  equivalent. 
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2.2  SEARCH  FOR  OBJECTS  ARRANGED  IN  A  PATTERN 


Now,  consider  the  construction  of  event  spaces  in  the  search  for  objects  placed  in  predefined 
patterns.  When  objects  are  placed  in  these  coordinated  positions,  they  no  longer  exhibit  the 
property  of  mutual  independence.  Indeed,  the  more  distinctive  the  pattern  is,  the  more 
interdependent  the  object  placements  become. 

Continue  with  the  paradigm  of  placing  the  set  M  of  Nm  known  objects  in  the  search  space. 
Now,  however,  the  general  form  of  the  joint  probability  density  function  is  expressed  as  a 
product  of  conditional  densities,  using  Bayes  rule: 


f(x,,...,xN  )  =  f,(x,)  nfj(Xj  |X,,...,XH,). 

j=2 


(39) 


Here  f ,  ( x , )  denotes  an  unconditional  initial  object  placement.  This  form  for  the  joint  density 
highlights  dependency  on  objects  previously  placed  in  the  pattern.  When  a  sequential  placement 
depends  only  on  the  previous  placement,  the  joint  density  takes  the  form  of 


f(x, 


xNM)  =  fI(xI)rifj(xJ|xJ.„). 

j=2 


(40) 


Such  placements  often  occur  when  objects  are  deployed  in  sequence  from  a  moving  platform, 
such  as  a  ship  or  aircraft. 

In  the  joint  probability  density  formulations  of  equations  (39)  and  (40),  the  object  ordering 
is  specific  and  indicative  of  the  applied  deployment  strategy.  Alternately,  placement  may  be 
structured  about  a  common  fixed  reference  point,  as  illustrated  in  figure  5.  The  reference 
position  may  constitute  a  general  drop  location  in  a  compound  experiment  in  which  the  structure 
is  placed  in  the  search  space  as  a  unit.  For  this  case,  there  are  2Nm  +  2  degrees  of  freedom 
realized  in  the  field  and  object  placement — two  degrees  of  freedom  for  the  general  drop  and  two 
for  each  of  the  Nm  objects  relative  to  the  unit  center.  The  joint  object  placement  density  is 
obtained  through  marginalization  over  the  field  placement  variable: 


f(x„...,xNM)  =  |s  f(x„...,xNM  |  xref)fref(xref)dxref, 


=  fs 


nfj(Xjixref) 

j=1 


fref(Xref)dXref 


(41) 


With  this  modeling  paradigm,  object  placement  is  conditionally  independent  within  the  structure 
but  reference  point  dependent  over  the  search  space. 
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a.  Sequential  Placement 


b.  Reference  Point  Placement 


Figure  5.  Object  Placement  Conditioning  Variables 


The  properties  of  this  object  placement  formulation  are  now  examined  in  more  detail.  The 
patterns  depicted  in  figure  5  are  characterized  by  a  separation  between  objects.  Assume  that 

object  placement  is  restricted  to  occur  over  a  set  of  regions  Q  =  jo  ,  ,  such  that 


x,  e  Qj,  Xj  g  Q|,  j  =  i  *  j .  (42) 

The  set  Q  constitutes  the  disjoint  regions  in  R2  over  which  the  reference  conditional  probability 
density  functions  fj(xj  |  xrct)  have  support;  that  is,  fj(xj  |  xrcf)  =  0  for  x )  g  Q  i . 


Assume  that  the  reference  location  xrcf  is  given.  Let  the  regions 

Rj  =  Qj  n  [J  Ai  >  (43) 

ielp 

for  j  =  1,...,Nm,  denote  the  intersection  of  the  conditional  support  region  with  the  partial  search 
path.  Then,  the  reference  conditional  clearance  probability  becomes 


Pr(Ec  |xref)=  P  Pr(Xj  e  U  A,  nd)  =  d|  xref ) , 


H 

Nv 


ielp 


=  nJR  Pr(4=d  xj)fj(xj|xlef)dxj. 

J=1 


(44) 


Here,  the  detection  probability  is  considered  to  be  independent  of  the  reference  point  placement: 


Pr(4  =d|Xj,xref)  =  P  r(4  =d|xj). 


(45) 


The  clearance  probability  is  obtained  by  taking  the  expectation  of  equation  (44)  with  respect  to 
the  reference  placement  random  variable  xrcf  as 
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(46) 


Note  that  the  latter  form  of  equation  (46)  operates  over  the  reference  marginalized  probability 
density  function.  In  similar  fashion,  the  conditional  probability  of  detecting  at  least  one  of  the 
objects  becomes 


"•••  ]fref(Xref)dXref> 

=  JUA  Pr(<f|  =  d|  x)[  X  fj(x)  |  dx-C, 

»ip  '  u=l 


where  the  product  of  reference  conditional  complement  probabilities  is  expanded  into  a  linear 
combination  of  joint  probability  calculations  of  the  form  used  in  equation  (46).  The  parameter  C 
denotes  the  remaining  terms  of  this  finite  alternating  series  beyond  the  additive  component. 
Hence,  the  latter  form  is  constructed  to  denote  an  unsealed  mixture  minus  a  non-negative 
correction  term.  For  a  given  search  path,  the  reference  conditional  probability  density  functions 
fj  (xj  |  xrcf)  can  change  significantly  with  reference  location.  In  general,  the  component  integrals 
in  equations  (46)  and  (47)  are  distinct  for  each  object  mj  in  M,  and  functional  composition  cannot 
be  applied  as  in  equations  (32)  and  (35). 

Consider  as  an  example  a  structure  comprising  Gaussian  components  placed  about  the 
common  reference  point  xrcf.  Let  ^  =  xref  +  AXj  designate  the  intended  position  of  object  mj 

with  displacement  Axt  from  the  reference.  Assume  that  the  structural  alignment  of  the  group  of 
objects  does  not  change  as  the  reference  point  xrcf  is  moved  within  S.  Let  each  object’s 
conditional  placement  be  distributed  as  Xj  ~  N[Xj,I,ap]  and  assume  that  ap  «  | O  ( | ,  such  that 

any  probability  mass  outside  Qj  is  negligible.  Then,  for  j  =  1,...,Nm, 


(48) 
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and  the  joint  conditional  density  takes  the  form 


Nm 


f(x,,...,xNJ  Xref)=  nfJ(xJ  |xref). 


j=l 


1 


[2™l] 


7INTexp 


1  nm 


2ct"  j=i 


(49) 


Now,  let  the  reference  placement  density  also  be  Gaussian,  with  xref  ~  N[x0,I2a^ef  ] .  Following 
the  derivation  presented  in  appendix  B,  the  unconditional  joint  density  is  given  by 


f(x,,...,xN  )  =  ■ 


1 

( 

Pp 

expj 

f  i  f 

[2™p] 

|Nm 

vNM  +pIj 

[  2a;  l 

x  exp< - 


2o; 


v  +  Pp  j 


N„  „  V 

ZxJ-xJ  +  Zy,  -y, 

J=i  /  V  J=> 


(50) 


where  p,,  =  ap  /  aref  represents  the  ratio  of  individual  object  placement  and  group  placement 

standard  deviations.  Note  that  the  x  and  y  Cartesian  coordinates  associated  with  any  given  object 
location  are  uncorrelated.  However,  the  correlation  between  any  two  distinct  object  locations  in 
either  coordinate  direction  is  given  by  the  coefficient 


P(X|,,X.  )  = 


CTj,j; 

°j,CTh 


ref 


CTP  +  CTrer 


1  +  PP 


(51) 


Such  correlation  highlights  the  loss  of  independence  in  object  placement  for  this  jointly  Gaussian 
context. 


Observe  from  equation  (5 1 )  that  when  » 1  the  correlation  between  object  locations  is 

very  small.  This  happens  when  field  placement  is  very  specific  or  when  its  location  is  inferred 
via  the  fusion  of  identified  object  locations.  For  this  case,  the  independence  of  component 
placement  can  be  assumed  and 


fj(X|)  MojO  +  l/pj) 


rscxp 


2ctp 


[2rc]a 


Tcxp 


2  a 


(T^(x,-x„-Ax,r 

(xj-x„-A  xj 


(52) 


However,  the  object  locations  are  not  identically  distributed.  The  distinctiveness  of  the 
component  densities  is  governed  by  the  size  of  the  reference  point  displacements  relative  to  the 

drop  uncertainty.  Identical  distribution  of  components  can  be  assumed  only  when  lAxJ  «  ap  . 
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Generally,  the  structure  is  assumed  to  be  composed  of  distinct  objects.  Conversely,  as  p“  « 1 , 

the  field  placement  uncertainty  dominates  and  the  referenced  marginalized  components  become 
highly  correlated  but  approximately  identical  in  form.  In  such  a  case,  the  search  problem 
devolves  into  one  in  which  the  structure  itself  becomes  the  object  of  the  search  rather  than 
individual  objects  within  the  structure. 

These  two  limiting  conditions  serve  as  a  basis  for  search  evaluation  under  a  restricted  type 
of  placement.  Specifically,  the  likelihood  function  for  a  single  object  placement  is  developed, 
where  knowledge  of  field  structure  dictates  the  shape  of  the  single  placement  density.  Within 
this  likelihood  context,  the  capability  to  form  accurate  predictions  regarding  multiple  detection 
events  is  lost.  However,  a  performance  assessment  on  the  next  detection  probability  can  still  be 
attained  for  search  path  comparison.  To  do  this,  the  projection  of  the  reference  conditional 
probability  density  function  onto  a  single  object  placement  space  is  modeled  as  a  mixture 
density.  Hence,  the  single  object  placement  paradigm  reverts  to  an  independent,  identically 
distributed  sampling  context.  Each  component  of  the  mixture  density  corresponds  to  one  of  the 
plausible  object  locations  defined  for  the  field  structure: 

fm(x|  x«f )  =  Sfmj(xlxref»mj)Pr(mj) .  (53) 

H 


This  reference  conditional  density  function  has  support  over  general  coordinates  x  e  S.  The 
unconditional  density  is  produced  via  marginalization  over  the  field  placement  variable,  as  in 
equation  (4 1 ),  to  yield 


fm(*)=fs  Zfn,,(XlXref’mj)Pr(mj) 

H 


fref^ref)*1* 


ref  ' 


(54) 


In  the  mixture  model,  Pr(mj)  denotes  the  mixture  weight  corresponding  to  the  occurrence 
probability  for  object  mj.  For  equation  (53)  to  represent  a  proper  probability  density  function,  it 
is  required  that 

XPr(mJ)  =  l,  (55) 

j=i 

which  is  achieved  through  a  renormalization  of  the  individual  object  occurrence  probabilities. 

The  density  function  in  equation  (54)  can  be  used  to  compute  both  the  clearance  probability 
and/or  the  probability  of  object  detection  in  a  manner  analogous  to  equations  (46)  and  (47), 
respectively. 

The  detection  events  for  this  construct  correspond  to  an  unstructured  field  placement.  That 
is,  an  object  can  occur  at  any  location  covered  by  the  mixture  density.  However,  the  formulation 
does  allow  for  multiple  instantiations  occurring  within  the  same  reference  field  component.  This 
is  especially  relevant  when  identical  distributions  are  assumed  in  placing  Nm  objects  in  the 
search  space.  As  for  the  unstructured  placement  case,  when  equation  (35)  is  applied  to 
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determine  the  next  detection  probability,  the  same  functional  composition  effectively  reduces  the 
numerical  computation  from  the  SNm  x  5)Nv<  event  space  to  the  S  x  iD  subspace. 

Next,  the  performance  assessment  on  the  next  detection  probability  conditioned  on  a 
previously  identified  object  is  developed.  The  probability  question  now  becomes  finding  the 
likely  location  of  objects  given  that  one  of  the  objects  has  been  detected  and  identified.  Let  us 
presume  that  it  is  unknown  as  to  which  of  the  ordered  objects  is  the  identified  one.  Let  xl(i 
denote  the  location  of  the  identified  object.  The  identified  object  location  is  related  to  the 
reference  point  location  via 

X,d  =  Xref  +  Ax,  +  Xei  =  *i  +  Xei  >  (56) 

where  xej  =  x-iq  ~  N[0,I2a^]  denotes  a  random  placement  error  (assuming  object  m,  as  the 
identified  object).  Assume  identical  error  distributions  fc(x)  for  each  component.  Then, 

Xj  =  Xref  +  AXj  =  X,d  +  AXj  -  AX,  -  Xei  .  <57) 

and  the  conditional  placement  density  function  for  object  irij  takes  the  form 

fmj(x|x,d,m1,rnJ)=  |s  fe(x-xid  -AXj  +  Ax,  +xcl)fe(xei)dxcj .  (58) 

These  density  functions  (for  the  various  m/s)  combine  in  a  mixture  form  as  in  equation  (53). 
Marginalizing  over  the  field  components  for  the  identified  (mi)  and  anticipated  (irij)  components 
yields  the  identified  object  placement  density  function: 


f«<xWs 


Ul  j=l 


fe(X-X,d>dx 


(59) 


Generally,  the  mixture  weights  are  not  dependent  on  object  location.  However,  if  the 
contribution  of  a  hypothesized  component  to  the  mixture  is  disallowed  to  enforce  object 
separation,  then  the  conditional  mixture  weights  must  be  renormalized  as 

i  Pr(m,)  .  .  .... 

Pr(m  m)  = - ^ ,  i*j  60) 

Jl  1  1  -  Pr(mj ) 

to  guarantee  a  proper  probability  density  function. 


2.3  SEARCH  FOR  AN  UNKNOWN  NUMBER  OF  OBJECTS 

When  the  actual  number  of  objects  for  which  search  is  conducted  is  unknown,  the  search 
performance  evaluation  must  apply  scrutiny  over  any  object  placement  information  that  is  given. 
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In  the  absence  of  any  prior  knowledge  on  the  number  of  objects  or  the  placement  strategy,  it  is 
left  to  evaluate  the  search  effort  alone.  Anticipatory  knowledge  on  placement  number  may  come 
in  different  forms,  one  being  the  expected  value  of  the  number  of  objects  placed.  To  use  this 
knowledge  in  developing  the  joint  placement  probability,  object  placement  is  modeled  as  a 
spatial  Poisson  process. 

A  review  of  spatial  Poisson  processes  is  provided  in  appendix  A  (section  A.2),  and  an 
excellent  overview  of  Poisson  processes  is  given  in  reference  10.  The  use  of  Poisson  processes 
in  search  theory  is  well  studied  (reference  1 1);  they  are  used  extensively  in  modeling  false 
alarms.  In  summary,  probability  mass  is  assigned  to  every  possible  value  of  the  non-negative 
integer  nM  according  to  the  known  intensity  function  A.(x).  More  specifically,  for  any  given 
AcS,  the  probability  of  finding  N(A)  objects  is  given  by 

Pr(N(A)  =  nJ  =  e-X(A)^^.  (61) 

When  A,(x)  is  not  constant  over  the  space  S,  placement  occurs  as  a  non-homogeneous  spatial 
Poisson  process.  The  coloring  theorem  (see  appendix  A,  section  A.2)  allows  us  to  “paint”  the 
occurrences  of  objects  within  the  space  according  to  any  independent  criteria,  such  as  detection 
or  classification.  When  objects  are  colored  according  to  detection  and  non-detection,  the  result  is 
a  pair  of  independent  spatial  Poisson  processes. 

The  detection  probability  recursions  provided  in  equations  (26)  and  (27)  for  multiple  scans 
of  a  given  location  can  also  be  applied  to  the  spatial  Poisson  process.  At  a  location  xeS,  the 
probability  of  detecting  an  object  after  n  previously  unsuccessful  search  passes  yields  the 
detection  coloring 


y(x,n)  =  Pn(x)  •  P5(x,n  - 1) , 


(62) 


where  n  corresponds  to  the  nlh  passing  of  location  x  along  the  search  path.  The  detection 
intensity  function  for  the  entire  search  path  coverage  follows  as 


n 


Mx)  =  X(x)£y(x,  j)  =  X(x)  •  yn(x), 
j=i 


Yn(x)  =  lY(Xd) 


j=l 


(63) 


The  subscript  n  notation  is  used  loosely  here.  The  intended  meaning  is  that  for  any  x,  n(x)  is 
incremented  with  each  succeeding  search  pass  that  includes  x,  and  spawns  yet  another  set  of 
independent  spatial  Poisson  processes  over  the  repeated  search  region.  Over  multiple  passes,  the 
superposition  theorem  is  applied  to  aggregate  the  intensity  functions  and  Xn(x)  denotes  the 
aggregated  intensity.  For  the  given  search  path,  the  probability  of  finding  objects  within  the  grid 
cell  Gk.i  then  takes  the  form 


=  n. 


nM‘ 


(64) 
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As  cells  form  disjoint  sets  in  S,  the  independence  of  each  N(Gk,i)  allows  us  to  develop  the 
aggregate  number  of  detections  N|p  encountered  over  a  partial  search  as 


Pr(N,  =  nM  )=  e 


(»■,> 


Mp  V  'P 


n,,! 


where 


-i 

=  MUA,)* 

X  MGy) 

X  X  (n(k,l))  •  X.k , 

ielp 

(k,l)e  UA, 

ielp 

i  e  IP  (k,l)el, 

(65) 


(66) 


and 


A,kl  =  f  A,(x)dx .  (67) 

'  jGk.i 

Note  that  these  expressions  assume  a  grid  cell  resolution  that  is  fine  enough  to  allow  assumption 
of  constant  coloring  within  the  grid  cell. 

As  in  equation  (35),  equation  (65)  provides  the  basis  for  calculating  the  probability  of 
finding  at  least  one  object  along  the  partial  search  path  as 

Pr(Nlp  >l)=l-Pr(Nlp  =o)=l-e  X|p  .  (68) 

The  Poisson  process  model  assumes  that  the  state  variables  Xj  are  identically  distributed,  leading 
to  the  equivalent  probability  representations  of  equations  (38)  and  (64). 

Next,  the  probability  of  detecting  m  objects,  conditioned  on  a  known  number  of  objects  nM 
in  the  space  S,  is  examined.  Observe  that,  for  a  given  number  of  objects  N(S)  =  nM  in  S, 

Pr(N(S)  =  nM  )  =  e^(S)  (69) 

N(S)  can  be  decomposed  into  its  detected  and  undetected  components: 

N(S)  =  Nd(S)  +  N5(S)  =  Nd((J  A, )  +  N5(S) ,  (70) 

ieip 

where  it  is  recognized  that  detections  can  occur  only  within  the  partial  search  path.  Similarly, 
the  object  intensity  function  X(x)  is  decomposed  into  its  detection  and  non-detection  components 
as 
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X(x)  =  yn(x)A.(x)  +  (1  -  yn(\))X(x)  =  Xn(x)  +  X-(x) , 


(71) 


X5(x)  =  X(x)-?in(x). 

As  the  detection  and  non-detection  processes  are  independent,  it  follows  that  the  probability  of 
detecting  m  objects  within  the  partial  search  path,  given  nM  objects  in  S,  becomes 


Pr 


Nd(UA,)  =  m|N(S)  =  n, 


f 

f  > 

\ 

Pr 

Nd 

U  Aj 

=  m 

V 

l,e'p  ) 

) 

Pr(N3(S)=n„-m) 


Pr(N(S)  =  nM ) 


\k(kY  s)-M,b”‘ 

/ 

'e-MO  (MS))’”  ' 

m!  (nM-m)! 

/ 

nM! 

n»!  (KY(KW-K)T'm 

m!(nr1  -  m)!  (MS))"" 


nM^ 

(M  1 

m 

f  M  ) 

U(S)J 

l  MS)J 

(72) 


As  expected,  the  expression  in  equation  (72)  is  the  same  as  equation  (37),  with  the 
probability  p  represented  by  a  ratio  of  intensity  function  evaluations.  Equations  (36)  and  (66) 
highlight  another  interesting  property  of  the  intensity  function  X.(x).  Recognizing  the  rate  N  of 
object  placement  in  the  space  S  as  the  integral 

X(S)=  Js?i(x)dx  =  N,  (73) 


and  then  normalizing  X(x)  by  this  factor  yields  the  probability  density  function 


fm(x) 


M*) 

JsMx)dx 


Hence,  the  intensity  function  ratio  of  equation  (72)  takes  on  the  form  of  equation  (37)  with 


(74) 


X, 

*P 

MS) 


N'JuA,Yn(X)'fm(X)dX 

_ i*ip _ 

N 


X  ip^Mk.m  Pc,, 

i€lp(k,l)elj 


(75) 


In  general,  a  ^-decomposition  can  be  performed  on  any  intensity  function  to  yield  a  probability 
density  function  for  single  object  placement. 

The  application  of  the  spatial  Poisson  process  above  enables  the  modeling  of  identically 
distributed  object  placement  when  the  actual  number  of  objects  place  is  unknown.  It  is  also 
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consistent  with  the  formulation  when  the  exact  number  of  objects  is  known.  The  required  input 
is  in  the  expected  rate  of  occurrence  over  the  space  S.  As  a  consequence,  the  search  evaluation 
procedure  can  model  spatial  variations  in  object  placement  and  sensor  detection  performance. 

The  model  can  also  ingest  prior  knowledge  of  field  structure  when  the  restricted  context  of  single 
object  placement  is  applied.  This  can  be  done  for  the  case  of  identified  objects  or  for  the  general 
application  of  drop  reference  and  drop  preference  priors. 

To  construct  a  general  object  density  function  to  be  applied  in  the  evaluation  of  search,  let 
us  start  with  a  random  mine  placement  component  X.R(x), 


XR(x)  =  NR-fR(x),  (76) 

with  an  expected  number  of  randomly  placed  objects  of  N  R  .  This  intensity  function  presumes 

independent  and  identically  distributed  object  placements.  Let  ^m(x)  denote  the  intensity 
function  associated  with  a  single  identified  object  located  at  xip.  Then, 

^|D(X)  =  NID-fID(X)’  (77) 

where  the  form  of  fm(x)  is  that  of  equation  (59).  The  parameter  NID  specifies  the  expected 
number  of  undetected  objects  remaining  within  the  field  structure.  As  noted  in  section  2.2, 
independence  can  be  assumed  if  the  positional  uncertainty  of  Xid  is  small  relative  to  the  field 
structure.  That  is,  the  next  detection  event  could  occur  anywhere  in  which  I'id(x)  has  support. 

The  identical  distribution  assumption  is  valid  as  long  as  detection  events  are  restricted  to  the  next 
detection.  If  additional  objects  affiliated  with  this  structure  become  identified,  then  Xid(x)  and 
fio(x)  would  need  to  be  updated  before  they  can  be  applied  again. 

Consider  the  field  prior  A.Mp(x)  associated  with  object  placement  in  the  space  S  given  a 
disposition  to  place  objects  in  the  designated  structure.  This  intensity  function  has  the  similar  X- 
decomposition 

^MF(x)  =  NMF-fMF(x).  (78) 

For  this  case,  the  parameter  NMF  can  take  on  several  meanings  and  its  value  must  be  set 
according  to  the  intended  assessment.  Let  fMF(x)  take  on  the  form  of  equation  (54).  As 
discussed  in  section  2.2,  while  the  reference  conditional  density  of  equation  (53)  is  conditionally 
independent,  a  large  drop  preference  region  will  induce  a  high  degree  of  correlation  between 
objects  that  are  placed  first  in  the  structure  and  then  mapped  to  S.  Thus,  independence  in 
placement  does  not  hold.  However,  independence  does  hold  in  the  placement  of  the  structure 
within  S.  Consequently,  NMF  can  be  set  to  indicate  the  expected  number  of  field  placements. 

The  actual  number  of  field  placements  may  or  may  not  be  known.  If  known,  then  the  density 
fMi  (x)  is  applied  directly  in  the  search  evaluation  to  assess  the  probability  of  finding  the 
structure.  If  unknown,  then  the  spatial  Poisson  process  of  equation  (78)  is  used.  Furthermore, 
field  presence  assessments  can  be  conducted  separately  from  next  object  detection  probability 
evaluations.  As  every  field  instantiation  must  contain  at  least  one  object,  limiting  values  for 
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N  MF  can  serve  as  a  bound  on  the  next  detection  probability.  An  approximate  measure  of 
performance  can  be  achieved  using  the  product 

N„F  =  NF-NFM,  (79) 

where  N  F  denotes  the  expected  number  of  fields  and  N  FM  denotes  the  expected  number  of 
objects  within  the  field.  Care  must  be  applied  in  deciding  how  to  set  these  parameters.  The  size 
of  the  field  structure  relative  to  that  of  the  drop  preference  density  function  and  the  size  of  S 
should  be  considered.  Development  of  the  methodology  to  set  this  number  as  a  function  of  the 
field  prior  content  is  the  subject  of  future  investigations  and  is  not  discussed  further. 

Lastly,  an  aggregate  intensity  function  can  be  developed  as 

/,<x>  =  /,,(x)  +  /,ID(x)  +  /,M[(x), 

=  N„  ■  f,O0  +  N,„  ■  fID(x)  +  Nmf  •  fMF(i)  =  N  •  fm(x). 

This  object  intensity  function  summarizes  general  object  placement  conditions.  Here,  N 
denotes  the  total  expected  number  of  objects  in  S  regardless  of  source. 


2.4  ADDITIONAL  CLEARANCE  CONSIDERATIONS 

Clearance  implies  the  absence  of  undesirable  objects  within  the  space  S.  The  probability  of 
finding  Nm  known  objects  is  given  by  equation  (3 1)  or  equation  (32),  depending  on  the 
distribution  of  objects.  While  clearance  implies  the  neutralization  and/or  removal  of  objects 
from  the  space,  for  the  purpose  of  search  evaluation,  it  is  assumed  that  detection  implies 
neutralization.  That  is,  the  probability  of  detection  becomes  the  clearance  probability  as  once  an 
object  is  found  it  is  incumbent  on  the  searcher  to  determine  how  to  neutralize  it.  More  generally, 
other  factors  need  to  be  considered  for  neutralization  probability.  Equation  (46)  provides  the 
clearance  probability  when  the  field  structure  is  known.  This  calculation,  while  specific  in  form, 
is  often  difficult  to  calculate  in  practice  and  may  lose  relevance  when  precise  knowledge  is 
unavailable.  Equation  (68)  provides  the  probability  of  next  object  detection  for  the  general  case 
of  an  unknown  number  of  objects.  As  the  expected  number  of  objects  in  S  diminishes,  this 
calculated  value  decreases.  In  the  absence  of  new  detections  as  the  search  plan  is  executed,  a 
threshold  can  be  applied  to  indicate  a  degree  of  clearance.  Should  an  object  become  detected, 
the  calculation  must  be  reset  and  the  search  conducted  over  again. 

Clearance  probability  can  also  be  defined  relative  to  the  search  space  S.  To  accomplish  this, 
the  ^.-decomposition  of  the  aggregate  spatial  Poisson  process  is  used  and  the  placement 
conditional  probability  density  function  Pg,  as  defined  by  equation  (20),  is  operated  on.  Recall 
from  equations  (25)  and  (27)  that  non-detection  ( P^- )  is  defined  globally  over  the  space  using  a 

rectangular-grid  partitioning  of  S.  The  calculated  value  of  this  probability  is  conditioned  on 
object  placement  within  the  grid  cell  and  the  number  of  independent  searches  conducted  over  the 
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cell.  Define  E^  as  the  event  that  an  object  exists  within  cell  Gk.i  and  is  not  detected,  and  define 
E  j  as  the  event  that  an  object  exists  in  S  and  is  not  detected.  Rewriting  equation  (27), 


(81) 


i=l 


where  NSk  denotes  the  number  of  times  that  the  cell  is  interrogated.  The  probability  of  event 
Ej  occurring  is  calculated  through  a  marginalization  over  S,  as 


(82) 


A  single  object  clearance  criterion  is  then  constructed  as  a  threshold  on  the  probability  of  the  set 
complement  of  this  event: 


(83) 


Equation  (83)  provides  the  probability  that  if  there  is  an  object  placed  within  the  space  S  it 
would  be  detected.  The  density  function  Pg  provides  the  basis  for  calculating  a  weighted 
average  of  the  individual  cell  non-detection  probabilities  based  on  the  spatial  variability  of 
sensor  characterizations.  When  Pg  is  insufficiently  defined,  a  diffuse  prior  on  placement  can  be 
applied,  yielding  the  unweighted  average.  As  a  statistic  of  the  detection  event,  Pr(  Ed ) 

represents  an  evaluation  of  the  search  effort  applied  over  S,  rather  than  a  measure  of  the 
likelihood  of  finding  an  object  within  S.  Greater  scrutiny  on  the  clearance  condition  can  be 
achieved  by  thresholding  the  criterion  at  the  grid  cell  level  in  lieu  of  the  aggregate  statistic.  For 
this  case,  if  the  threshold  probability  is  not  attained  anywhere  within  the  grid,  then  it  is  not 
attained  over  the  space.  This  is  equivalent  to  thresholding  the  probability  test  statistic: 


(84) 
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3.  SEARCH  PATH  EVALUATION 


Drawing  heavily  on  the  development  of  probability  event  spaces  presented  in  the  preceding 
section,  this  section  describes  the  design  of  a  grid-based  search  path  evaluation  capability.  The 
general  requirements  for  this  capability  are  fairly  simple — the  methods  must  be  general  enough 
to  support  evaluation  of  a  variety  of  alternate  search  path  generation  techniques  yet  specific 
enough  to  provide  realistic  assessments  for  scenarios  that  exhibit  a  high  degree  of  contextual 
complexity.  For  the  mine-hunting  application,  this  includes  the  adoption  of  spatially  varied  mine 
and  background  clutter  densities  and  environmentally  dependent  sensor  performance 
characterizations.  The  evaluation  capability  must  be  applicable  when  the  number  of  mines  in  the 
search  region  is  either  known  or  unknown.  In  addition,  the  capability  must  support  search 
evaluation  in  situations  with  uncertain  descriptions  of  mine  field  placement  structure. 


3.1  GENERAL  SETUP  OF  THE  EVALUATION  SPACE 

Let  ScR2  represent  a  physical  region  in  geographic  coordinates  and  over  which  the  search 
for  mines  is  to  be  conducted.  A  grid  is  developed  over  these  spatial  coordinates  to  a  resolution 
that  accommodates  the  natural  variability  in  mine  placement  strategies  that  occur  due  to 
dependencies  on  bathymetry,  bottom  composition,  and  geographic  context  (such  as  the  proximity 
to  harbors  and  shipping  lanes).  For  simplicity  of  illustration,  a  square  grid  space  is  constructed 
with  a  rectangular  tessellation,  as  illustrated  in  figure  6.  However,  the  grid  need  not  take  a 
rectangular  shape  nor  must  the  space  be  uniformly  sampled.  This  simpler  realization  places 
emphasis  on  development  of  the  underlying  likelihood  structure  over  which  the  probabilistic 
inferencing  is  conducted.  Probability  calculations  are  performed  over  the  grid 
G  =  { x  =  (x,y)  €  S|  x  =  ±(k  -1/2)AX,  y  =  ±(m-l/2)AY,  (k,m)  =  l,...,Nmax}  defining  the  set  of 

geographic  coordinates  in  a  square  of  side  length  2NmaxAX  and  equal  spacing  AX  =  AY  centered 
about  (x,  y)  =  (0,0) .  Hence,  4Nmax2  grid  points  are  used  to  construct  the  space.  Probabilities  are 

updated  sequentially  at  a  time  sampling  that  is  specified  to  occur  at  regular  intervals 
T  =  { tj  |tj  =  t0  +  iAT,i  =  0,l,...,Imax },  where  time  to  indicates  the  start  time  of  the  search.  The 

sampling  interval  AT  is  constrained  to  be  sufficiently  small  relative  to  vehicular  dynamics  such 
that  the  search  trajectory  over  this  interval  is  approximately  linear.  However,  the  constraint  ot  a 
small  and  uniform  AT  can  be  relaxed  in  alternate  realizations  without  altering  the  likelihood 
structure. 

In  figure  6,  the  square  search  region  appears  on  the  left,  with  mine  likelihood  contours  and 
sample  objects  (mines  and  clutter)  overlaid  on  it.  The  mine  locations  are  represented  by  circles 
and  the  clutter  locations  by  plus  signs.  A  segment  of  a  random  search  path  is  also  illustrated  in 
the  figure.  The  end  of  the  search  path  that  contains  a  circle  and  box  indicates  the  sensor  location 
at  a  current  time.  The  box  defines  the  region  within  the  space  used  to  construct  a  local  sensor 
referenced  map  of  the  following  form: 

L  =  {x  |  x  eG,  xs  -dx  <  x  <  xs  +dx,  ys  -dy  <  y  <  ys  +dy}-  (85) 
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Figure  6.  Evaluation  Grid  Development 


This  local  map  centers  at  the  current  (time  =  tj)  sensor  coordinate  reference  Xs(ti),  as  is  illustrated 
in  the  lower  right  portion  of  figure  6.  Here,  the  uniform  grid  spacing  and  the  square  cell 
construction  of  G  becomes  evident.  Displacements  dx  and  dy  are  selected  to  ensure  that,  for  the 
given  update  interval  AT,  all  points  covered  by  the  sensor  scan  along  this  trajectory  segment  are 
included  in  L.  A  linear  transformation  of  coordinates  in  the  local  map  at  time  tj  of  the  form 


xMi(tl)  =  T(0)[x-xs(ti)l  T(0(t,)) 


0056(11)  -sin  0(f) 
sin  0(f)  0056(11) 


(86) 


is  applied  to  align  coordinates  with  the  vehicle  orientation  0(tj).  This  is  illustrated  in  the  upper 
right  portion  of  figure  6.  The  circles  are  of  radius  Rd,  which  is  indicative  of  the  detection  region 
covered  by  the  sensor.  Also  indicated  is  the  change  of  coverage  as  the  sensor  advances  from 
time  tj  to  tj+i  along  its  prescribed  search  path.  The  apparent  slither  at  the  top  of  the  circle  is 
defined  as  the  evaluation  region. 
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Within  the  evaluation  process,  multiple  maps  are  developed  over  these  global  geo-spatial 
coordinates,  including  a  mine  density  map,  a  clutter  density  map,  and  a  geographically 
referenced  sensor  detection  map.  Each  of  these  maps  associates  probability  mass  to  the  grid 
elements,  and  all  of  the  maps  conform  to  the  same  grid  specification  G.  The  mine  density  map 
forms  the  cornerstone  on  which  search  effectiveness  evaluation  is  conducted.  As  discussed  in 
section  2.3,  the  mine  density  map  is  decomposed  into  a  probability  density  function  scaled  by  the 
expected  number  of  mines.  Globally  referenced  detection  modeling  is  used  to  isolate 
geographically  referenced  sensor  performance  dependencies  and  to  maintain  the  bookkeeping 
associated  with  multiple  pass  search  paths.  However,  all  probability  calculations  are  made  over 
the  local  sensor  referenced  map  (over  grid  L). 

The  search  region  size,  grid  resolution,  and  update  interval  can  be  specified  directly  or  can 
be  derived  as  a  function  of  three  control  variables — pA,  Nmax,  and  Nadv  In  the  latter  case,  as 
above,  let  Nmax  specify  the  number  of  grid  points  in  a  half-length  along  the  side  of  the  search 
region.  Then,  setting  each  grid  point  at  the  center  of  each  identical  grid  cell  of  area  Ax*Ay,  the 
size  of  search  region  S  becomes  ao  =  |S|  =  4Nmax2*Ax*Ay.  Let  pA  specify  the  ratio  of  sensor 
detection  radius  Rd  to  the  distance  advanced  in  a  search  update  interval  AT.  Then,  pA  =  RD/Dadv 
«  Rd/SsAT,  where  DadV  *  SsAT  for  a  search  platform  traveling  at  constant  velocity  (it  is  assumed 
that  the  update  interval  AT  is  small  enough  for  such  an  approximation  to  hold).  Let  NadV  be  a 
numerical  resolution  parameter  specifying  the  number  of  grid  points  in  the  search  advance  slither 
in  the  direction  of  forward  motion.  Thus,  the  temporal  update  interval  is  given  as  AT  = 

[l/pApD]f,  where  pD  =  Ss/Rd  and  the  plus  brackets  imply  rounding  up  to  the  nearest  integer 
value.  Similarly,  the  grid  spacing  resolution  in  either  Cartesian  coordinate  direction  is  calculated 
as  Ax  =  Ay  =  [Dadv/  Nadv]"  =  [Rd  /  PaN^v]",  where  the  minus  brackets  imply  rounding  down  to 
the  nearest  integer. 

This  grid  specification  highlights  interdependence  between  resolution  and  the  sequential 
update  interval  determined  as  a  function  of  the  sensor  characteristics  (Rd),  search  rate  (Ss  and 
AT),  and  the  mine  field  characteristics  (Ax,  Ay,  and  ao)  that  are  assumed.  Otherwise,  the  size  of 
the  search  region  may  be  specified  and  the  resulting  number  of  grid  cells  necessary  to  attain 
accuracy  determined.  That  is,  either  Nmax  or  ao  can  be  made  the  independent  variable. 

3.2  NUMERICAL  CALCULATION  OF  SEARCH  PROBABILITIES 

Recall  that  the  detection  event  decomposition  described  in  section  2. 1  requires  that,  for  each 
grid  cell,  a  corresponding  value  for  the  probability  of  detection  conditioned  on  arbitrary  mine 
placement  within  that  cell  be  calculated  as 


(87) 


where  |Gk,i|  is  the  area  of  the  grid  cell,  as  given  by  equation  (11).  The  form  of  equation  (87) 
assumes  constant  mine  placement  likelihood  over  the  grid  cell,  and  this  assumption  is  maintained 
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henceforth.  As  noted,  the  resulting  detection  probability  represents  a  spatially  averaged  value  for 
detection  probability.  When  the  detection  likelihood  within  the  grid  cell  is  also  approximately 
constant  (at  value  Pd),  equation  (87)  simplifies  to 


(88) 


as  expected.  Thus,  under  an  assumption  of  constant  detection  likelihood  within  the  grid  cell,  the 
grid  cell  conditional  detection  probability  can  be  assimilated  directly  from  a  suitable  sensor 
performance  prediction.  Note  that  this  expression  does  not  depend  explicitly  on  sensor  position. 
Rather,  sensor  proximity  was  implied  in  the  designation  of  search  regions  and  their 
correspondence  to  scans.  Digressing  from  this  restrictive  association,  the  paragraphs  below 
introduce  the  circular  sensor  evaluation  region  that  is  depicted  in  figure  6. 

Let  the  sensor  detection  capability  be  characterized  by  the  definite-range  law  (reference  3) 
with  a  fixed  and  known  detection  radius.  That  is,  if  the  target  appears  anywhere  within  a  circle 
of  radius  Rd  centered  at  the  sensor,  the  target  will  be  detected  by  that  sensor  with  probability  Pd. 
Consecutive  circular  regions  are  applied  for  each  temporal  update,  and  the  set  difference  as  the 
circular  region  at  time  tj.i  sweeps  to  that  at  time  t;  is  determined.  The  set  difference  constitutes 
what  is  referred  to  as  the  evaluation  region,  and  it  is  developed  assuming  linear  motion  between 
updates.  Hence,  its  shape  remains  constant  even  as  the  search  platform  maneuvers  about  in 
executing  its  search  plan.  From  this  region,  a  probability  change  AP  is  calculated  for  each  time 
interval  and  mapped  to  the  event  spaces  under  consideration. 

A  search  path  constitutes  one  contiguous  swath  centered  on  the  vehicular  trajectory 
established  by  the  platform  search  plan.  The  search  region  is  segmented  along  the  trajectory  at 
sampling  interval  AT  =  tj  -  tj_i.  This  is  illustrated  in  figure  7  for  an  arbitrary  search  path.  For 
non-intersecting  paths,  the  segmented  search  path  constitutes  a  collection  of  spatial  subsets  that 
are  both  exclusive  and  exhaustive.  This  makes  the  collection  well  suited  to  set  theoretic 
probability  calculations.  Note  that  the  segments  need  not  relate  directly  to  specific  sensor  scans. 
All  that  is  required  is  that  values  assigned  to  the  geographically  referenced  Pd  map  be 
representative  of  those  obtained  using  anticipated  sensors  traveling  along  the  search  trajectory. 

If  need  be,  the  Pd  map  can  be  developed  more  specifically  as  a  function  of  source-receiver 
geometry.  For  this  case,  proximity  to  the  search  trajectory  along  the  segment  must  be 
considered.  Given  a  complete  search  path,  the  Pd  map  can  be  constructed  over  the  entire 
geographic  space  prior  to  evaluation.  However,  altering  the  path  may  require  that  the  map  be 
updated  as  well. 
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Figure  7.  Example  of  Search  Path  Decomposition 


For  purposes  of  this  report,  the  distribution  of  detection  likelihood  within  the  search  path  is 
of  concern  only  for  the  case  of  intersecting  search  paths  or  when  the  placement  density  varies 
significantly  within  the  search  path  segment.  Otherwise,  applying  an  averaged  Pd  to  the  path 
yields  the  same  probability  calculation.  That  is,  if  mine  density  is  approximately  constant  over  a 
region  Grc  comprising  Nrc  grid  cells  in  the  set  Rc  as  in 

G»t  =  U°u  •  <89> 

<k.l)eRc 


then 

f  Pr(d  =  d  |Gk  ,)dx  f  Pr(d  =  d  |GRc)dx 

p  =— l—  y p  =  y  - .  (90) 

rDRc  N  Z-rDk.,  Z  NG  Ig 

'  Rc  (k,l)eRc  (k,l)eRc  1NRc|Uk,l|  I^Rcl 

The  current  realization  requires  the  development  of  the  global  Pd  map  to  facilitate  comparison  of 
alternate  search  paths.  However,  this  design  requirement  can  be  relaxed  should  detection 
probability  need  recalculation  on  a  more  frequent  basis.  For  such  a  case,  the  local  Pd  map  would 
be  developed  and  applied  within  each  update  cycle  in  lieu  of  a  global  articulation  prior  to 
evaluation. 


For  the  cases  of  intersecting  search  paths  and/or  multiple  search  platforms,  non-detection 
probability  must  be  maintained  at  a  global  scale  to  preserve  search  history.  Specifically,  a  non¬ 
detection  Pp  map  is  constructed  over  the  same  grid  G  and  initialized  to  unity  prior  to  the  start  of 

the  search.  Then,  for  each  time  that  a  grid  cell  Gk,i  is  covered  by  a  search  segment,  the  non¬ 
detection  map  is  updated  at  the  cell  index  k,l  according  to  the  recursion 


=  Pok,0-PDkl) 


(91) 
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for  all  the  points  in  the  evaluation  region.  Maintaining  this  quantity  separately  leaves  the 
detection  PD  map  unchanged.  The  multiple  pass  detection  probability  is  thus  provided  by  the 
PDP5  product  with  initial  value  PD. 


The  probability  change  AP,  over  a  time  interval  AT,  is  calculated  by  aggregating  the  joint 
probability  for  each  point  in  the  evaluation  region.  Let  Im  represent  the  set  of  indices  covered  by 
the  sensor  and  its  evaluation  region  at  time  tj.j.  The  set  I,  is  constructed  to  include  all  cell  indices 
covered  by  the  union  of  the  circular  sensor  coverage  region  at  time  tj  and  the  tubular  swath 
connecting  it  to  the  previous  sensor  coverage  region.  Then,  the  set  difference 

I Aj  =  {(M)  |(k,l)  e  Ij,(k,l)  I,  , }  (92) 


provides  the  indices  of  all  grid  points  that  lie  within  the  evaluation  region  developed  at  time  tj. 
In  this  scheme,  fringe  points  may  occur  arbitrarily  within  any  particular  update.  However,  the 
fringe  points  are  aggregated  only  once  per  distinct  search  path,  and  so  the  integrity  of  the 
calculation  is  preserved  over  paths  of  arbitrary  length  to  within  the  fringe  accuracy.  The 
probability  change  over  the  time  interval  follows  as 


AP,=  X  Pd  PB  -Pc 


(k,l)eIAl 


(93) 


The  quantity  P(,  represents  the  probability  density  function  provided  in  the  prior  knowledge 
specification  or  extracted  from  the  decomposition  of  the  mine  density  function  when  the  number 
of  mines  is  unknown.  Hence,  AP,  represents  the  probability  of  detecting  a  mine  within  the 
evaluation  region  conditioned  on  there  being  at  least  one  mine  within  the  search  space  S. 

For  the  case  of  an  unknown  number  of  mines,  the  mine  density  function  is  given  by 
^•m  =  N  m  "  Pq  ,  where  NM  denotes  the  expected  number  of  mines  in  the  search  region  S.  A 

value  for  NM  must  be  provided  in  the  prior  specification.  Then,  the  expected  number  of  mine 
detections  in  the  evaluation  region  becomes  NM  =  NM  •  AP, ,  and  the  probability  of  not  finding 
any  mines  in  the  evaluation  region  becomes 

AQss,  =exp[-NM  AP,]  .  (94) 


The  aggregate  non-detection  and  detection  probabilities  over  the  partial  search  path,  up  to  and 
including  segment  i,  become 


Qss,  —  FI  ^Qss.  —  ^Qss,Qss, 

j-0 

Pss,  :  1  —  Qss,  ’ 


(95) 
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respectively.  Analogous  quantities  are  developed  for  the  case  of  the  clutter  density  map 
Xc  =  Nc  •  PGc ,  where  the  value  of  Nc  — the  expected  number  of  clutter  objects  in  S — is  also 
assumed  to  be  provided  in  the  prior  knowledge  specification. 

For  the  case  when  Nm — the  actual  number  of  mines  in  S — is  assumed  known,  and  each 
mine  is  assumed  to  be  identically  distributed  with  density  function  Pq,  the  probability  calculation 
takes  the  forms  developed  in  section  2.1.  In  such  a  case,  the  clearance  probability  for  finding  all 
Nm  mines  in  the  partial  search  path  is  calculated  as 


P  = 

rCP, 


i  An 

2>p, 

Vj=° 


(96) 


Similarly,  the  probability  of  finding  at  least  one  of  the  Nm  mines  along  the  partial  search  path  is 
calculated  as 


Pmp,  =1- 


-5>p, 


\  Nm 


j=0 


(97) 


3.3  EXTENSIONS  FOR  MINE  FIELD  STRUCTURE 

Any  prior  knowledge  of  mine  field  structure  is  integrated  into  this  search  evaluation 
capability  in  the  generation  of  either  the  underlying  mine  density  map  XM  or  its  component 
probability  density  Pg-  The  fundamental  construct  is  the  notion  of  a  local  mine  field  map  Pmf, 
developed  over  a  scale  commensurate  with  the  prior  knowledge  of  field  structure.  This  field  map 
can  be  employed  in  a  variety  of  ways  in  calculating  event  probabilities.  The  exploitation  of  the 
field  maps  is  the  subject  of  ongoing  investigation.  The  following  representations  apply  to  the 
restricted  set  of  probability  calculations  whereby  identical  distribution  in  single  mine  placement 
can  be  assumed. 

The  mine  location  density  function  XM  is  constructed  as  a  superposition  of  an  unstructured 
mine  placement  density  and  the  density  function  associated  with  each  candidate  mine  field 
structure: 


-  ^-R  + 

j=> 


MF, 


(98) 


where  it  is  recalled  that  NM  denotes  the  expected  total  number  of  mines  in  the  search  region  S. 
Each  component  of  ^m  is  similarly  decomposed  as 


^r=Nr 


•P, 


GR  ’ 


^  MF;  —  N  Ut  ■  P< 


MF, 


1  GMF,  ’ 


j  =  1,...,N 


MF  ’ 


(99) 
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where  NR  is  the  expected  number  of  randomly  placed  mines  and  NMF  is  the  expected  number 

of  mines  placed  using  mine  field  structure  i.  Note  that  NMF  may  be  more  or  less  than  the 

number  of  mines  necessary  to  totally  populate  the  structure  N  MF  .  This  can  happen  when  the 

structure  population  is  unknown  or  when  multiple  realizations  of  the  field  can  be  expected.  For 
some  event  criteria,  the  value  may  signify  only  the  number  of  mine  fields  in  the  space  S. 

Substituting  equations  (99)  into  equation  (98)  and  rearranging  the  terms  allows  the 
probability  density  Pg  to  be  formulated  as  a  mixture  density  of  the  form 


P  =  W 

rG  VVR 


lGR 


Mt- 

2!  Wmf, 


1  GMF, 


j=l 


(100) 


where  mixture  weights  are  calculated  as  WR  =  Nr/Nm  and  WMF  =  Nmf  /Nm  .  In  this  way,  Pg 

incorporates  prior  knowledge  of  mine  field  structure  at  a  degree  proportional  to  its  anticipated 
occurrence.  The  field  structure  is  developed  in  the  form  of  a  conditional  probability  density 
function  as  described  in  section  2.2.  Specifically,  a  drop  referenced  mine  field  map  with  a 
central  drop  reference  point  is  constructed  in  local  mine  field  coordinates.  With  application  of 
the  drop  independence  assumption,  this  function  is  convolved  with  a  drop  preference  map 
developed  in  global  coordinates  to  obtain  the  resulting  mixture  component  PGMF  .  A  separate 

drop  referenced  mine  field  map  is  constructed  for  each  candidate  structure.  Drop  preference 
maps  may  be  specific  to  mine  field  structure,  depending  on  its  composition  and  any  constraints 
imposed  on  mine  types  included.  The  grid  resolution  for  the  drop  reference  and  drop  preference 
maps  is  identical  to  that  of  the  global  evaluation  grid  G.  However,  the  size  of  the  map  is  specific 
to  the  structure.  For  illustration,  a  construction  example  is  provided  below. 

Consider  an  example  in  which  the  parameterization  of  prior  knowledge  on  field  structure 
consists  of  a  single  mine  field  type  with  a  number  of  mine  lines  NL  =  2.  Each  line  consists  of  a 
fixed  number  of  mines  with  known  spacing  Dm  between  them.  The  spacing  Dl  between  mine 
lines  is  assumed  to  be  fixed  and  known.  For  this  example,  two  random  placement  influences  are 
introduced.  First,  a  Gaussian  mine  drop  error  of  the  form  xc  ~  N[0, 12CT2]  is  applied  for  each  mine 
in  the  structure.  In  this  example,  the  drop  error  variance  a2  is  circularly  symmetric,  with  a  = 
0.13*Dm-  The  second  random  influence  is  introduced  in  an  ancillary  parameter  that  is 
independent  of  the  space  S.  Define  the  angle  0  to  designate  the  orientation  of  the  mine  field 
(assuming  the  rotation  occurs  about  the  reference  point).  Let  Iq  =  [0O  -  A0/2,  0()  +  A0/2] 
represent  an  angular  interval  A0  about  0O  from  which  the  orientation  angle  0  is  drawn. 

Nominally,  all  ancillary  random  variables  must  be  marginalized  out  in  the  construction  of  the 
mine  field  map.  When  available,  statistical  parameterizations  of  field  structure  are  applied 
directly  in  the  marginalization  process.  Otherwise,  maximum  entropy  methods  (reference  1 2) 
are  applied  to  develop  plausible  mine  locations.  These  entropic  priors  require  a  minimal  amount 
of  modeling  assumptions  in  their  development.  An  example  is  illustrated  in  figure  8a.  The 
functional  form  used  in  developing  the  entropic  drop  referenced  map  shown  in  figure  8a  is 


fm(X)  =  fMF(XlXref)*fref(XrerK 


(101) 
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where 


fMF(xix 


ref 
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( 

max' 

exp 

u. 

_ 1 

V 

dp(x,xJ;0)  +  (pde(x>xj;0)) 
2ct2 


2\1 


(102) 


In  equation  (102),  Kmf  is  a  normalization  constant  to  ensure  that  equation  (101)  serves  as  a 
proper  probability  density  function.  Imf  indicates  the  indices  of  mines  within  the  structure  and  dp 
and  de  denote  radial  and  angular  distances  of  the  form 

dp(x,,x2;xref)  =  |p2  —  p,|,  (103) 

de(x,,x2;xref,A0)  =  max(max(02  - A0/2 -0,,O),max(02  +  A0/2-0,,O)),  (104) 

where  p,  =  ^/(x,  - xref )2  +  (y,  -yref)2  and  p2  =  ^/(x2  - xref )2  +  (y2  - yref )2  denote  the 
distances  from  the  reference  point  (supplied  as  an  ancillary  parameter)  to  the  points  X|,  X2,  and 
0,  =  tan  1  (y ,  -  yref/x,  -  xref )  and  02  =  tan  '(y2  -  yref/x2  -xref )  denote  the  angles  between  Xi, 

Xj  and  the  coordinate  axis.  The  radial  distance  of  equation  (103)  measures  the  absolute 
difference  between  the  two  calculated  radials  and  the  angular  distance  of  equation  ( 1 04) 
measures  the  difference  between  the  two  calculated  angles  relative  to  the  size  of  the  interval  A0. 
The  latter  comprises  an  angular  interval  distance.  Since  the  drop  reference  map  is  centered  about 
the  reference  point  location,  the  reference  location  used  in  equations  (103)  and  (104)  is  at  the 
coordinate  origin  of  the  mine  field  map,  xrcf  =  O2. 


a.  Drop  Reference  Map  b.  Drop  Preference  Map  c.  Drop  Referenced  Prior 


Figure  8.  Drop  Referenced  Prior  Development 
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An  example  of  a  drop  preference  map  is  shown  in  figure  8b,  and  the  composite  drop 
referenced  prior  is  generated  via  a  numerical  convolution  of  the  drop  reference  and  drop 
preference  maps.  The  result  is  shown  in  figure  8c.  Figure  9  illustrates  the  difference  in  scale 
between  the  search  space  and  the  mine  field  drop  reference  map.  The  size  of  the  mine  field  map 
is  indicated  by  the  blue  box  in  the  figure.  Overlaid  are  a  replica  of  the  drop  referenced  prior  and 
a  sample  instantiation  of  a  mine  field  consistent  with  the  assumptions  on  the  prior.  Background 
clutter  is  depicted  as  the  blue  and  red  plus  (+)  symbols,  with  red  indicating  a  mine-like 
appearance.  The  black  star  (*)  in  figure  8b  represents  a  numerically  drawn  sample  from  the  drop 
preference  map.  The  sample  mine  field  instantiation  is  centered  at  this  location. 


Figure  9.  Example  of  Drop  Referenced  Prior 


A  distinct  set  of  maps  are  generated  when  mines  are  known  to  be  present  in  the  search 
environment  a  priori.  This  is  so  because  of  the  extended  mixture  hypotheses  required  to  form 
the  map  as  shown  in  equation  (59).  Figure  10  illustrates  a  sample  set  of  the  maps.  The  mine 
location  map  of  figure  10a  represents  the  uncertainty  in  the  identified  mine  location.  This  can  be 
merely  a  numerical  evaluation  of  the  drop  error  uncertainty  or  can  incorporate  additional  factors 
such  as  registration  uncertainty  due  to  navigational  errors.  In  the  current  realization,  the  scale  of 
this  map  is  set  to  that  of  the  global  coordinate  grid  G.  In  practice,  this  function  has  small  support 
over  S  and,  therefore,  requires  a  grid  that  is  a  fraction  of  its  current  size.  The  global  scale  does 
accommodate  multiple  identified  mines  within  the  current  processing  methodology,  provided 
that  they  do  not  originate  from  the  same  field  structure.  The  plausible  location  map  in  figure  10b 
indicates  where  mines  might  be  found  given  an  identified  mine  location.  This  is  a  local  mine 
field  map  centered  on  the  identified  mine  coordinates.  The  current  realization  of  this  map 
applies  only  to  the  case  of  a  single  identified  mine.  Multiple  mine  identifications  that  may  have 
originated  from  the  same  field  structure  require  a  fusion  process  to  develop  a  posterior  location 
map  that  is  representative  of  the  multiple  observations.  The  map  developed  considers  each  mine 
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drop  location  represented  in  the  mine  field  structure  as  the  possible  position  of  the  identified 
mine.  Then,  for  each  of  these  points,  the  likelihood  of  other  mine  locations  within  the  mine  field 
map  is  evaluated  and  aggregated. 

For  purposes  of  illustration,  the  example  for  figure  10b  utilizes  an  entropic  parameterization 
in  the  map  formation  in  a  manner  similar  to  the  drop  reference  prior  example,  such  that 


—  fMp(X  I  XID)*fID(X.D)> 

and 


*Mf(X  I  X  ID  )  K-id 


max  < 

exp 

dp(x,xj;xi)  +  (pd9(x,xj;xi)) 
2cr 


(105) 


(106) 


where  denotes  the  normalization  factor.  Radial  and  angular  distance  metrics  are  developed 
as  in  equations  (103)  and  (104).  However,  the  reference  location  xrcf  now  becomes  the  mean 
location  for  the  hypothesized  identified  mine  component;  that  is,  xs .  The  convolution  of  the 
mine  location  and  plausible  location  maps  is  provided  in  the  field  structure  prior  illustrated  in 
figure  10c. 


a.  Mine  Location  Map  b.  Plausible  Location  Map  c.  Field  Structure  Prior 

Figure  10.  Identified  Mine  Prior  Development 


The  larger  scale  associated  with  the  field  structure  prior  is  illustrated  in  figure  1 1 .  The 
global  search  space  S  is  shown  on  the  right  and  a  blow-up  of  the  field  structure  prior  is  shown  on 
the  left.  The  identified  mine  location  was  set  to  the  first  instantiated  mine  of  a  sample  mine  field 
drawn  from  the  drop  reference  prior  of  the  previous  example.  The  sample  mine  field  is 
discernible  in  the  blow-up.  Observe  the  apparent  lack  of  likelihood  at  the  identified  mine 
position.  This  illustrates  the  point  that  the  mine  field  structure  is  modeled  to  preclude  mine 
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placements  in  close  proximity  to  each  other.  It  is  also  a  consequence  of  the  spacing  included  in 
the  prior  specification.  Note  in  equation  (100)  that  the  contribution  of  each  mine  field  map  is 
weighted  and  aggregated.  This  is  the  case  for  the  field  structure  prior  as  well.  Hence,  if  the 
expected  number  of  mines  N  Mh  is  on  the  order  of  the  other  priors,  the  field  structure  prior  can 

dominate  the  resulting  mine  density  map,  as  occurs  in  figure  1 1 .  The  example  depicts  the  prior 
considering  only  one  mine  field  type.  Multiple  types  can  be  aggregated  accordingly.  However, 
the  reference  for  any  composite  mine  field  structure  remains  at  the  identified  mine  location.  The 
determination  of  suitable  weighting  between  identified  mines  and  anticipated  mine  field  structure 
is  part  of  the  prior  knowledge  specification. 


Figure  1  I.  Example  of  Field  Structure  Prior  for  Identified  Mines 
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3.4  PROCESSING  DIAGRAMS 


The  following  diagrams  detail  the  sequence  of  processing  steps  in  the  computational 
implementation  of  the  search  evaluation  capability.  Figure  12  depicts  the  sequential  evaluation 
loop  that  is  realized  in  the  algorithm.  As  indicated,  results  are  provided  as  a  function  of  time 
along  the  search  path,  and  summary  statistics  are  generated  at  the  end  of  a  run.  The  two  essential 
processing  components  are  the  initialization  block  and  the  processing  cycle  conducted  for  each 
sequential  update. 


Figure  12.  Sequential  Evaluation  Loop  for  Search  Evaluation 


In  figure  13,  the  process  initialization  block  is  described  in  more  detail.  In  particular,  all 
prior  knowledge  of  field  structure  and  placement  is  applied  in  this  block.  The  result  of  this  prior 
knowledge  is  contained  in  the  form  of  the  mine  and  clutter  density  maps.  The  global 
representation  of  sensor  performance  is  developed  here  as  well.  Furthermore,  a  search  trajectory 
characterization  process  is  provided  within  this  block  to  synthesize  search  plans  for  basic 
motions.  This  component  employs  the  goal-oriented  curvilinear  motion  (GOCM)  model  that  is 
described  in  appendix  C.  In  practice,  most  search  plans  are  provided  as  direct  input  to  the 
evaluation  program  and,  in  such  cases,  this  block  reads  that  search  plan  into  memory. 

Figure  14  presents  a  detailed  description  of  the  update  cycle  itself.  Here,  the  search  path 
kinematic  trajectory  and  the  various  evaluation  maps  are  provided  as  input  and  operated  upon. 
The  diagram  illustrates  the  simplicity  of  the  modeling  construct.  Any  processing  component 
from  which  there  are  no  outgoing  flow  lines  indicates  the  availability  of  path-specific  results. 

The  branching  of  processing  components  illustrates  opportunities  for  calculations  to  be 
performed  in  parallel.  Following  completion  of  the  evaluation,  summary  results  are  produced  by 
either  extracting  the  path-specific  results  at  the  terminating  conditions,  or  by  operating  on  them 
in  some  ancillary  fashion. 
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Figure  13.  Initialization  Block  Diagram 
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Figure  14.  Update  Cycle  Diagram 


4.  ALGORITHM  PERFORMANCE 


This  section  presents  results  that  serve  to  validate  the  search  evaluation  methods  described 
in  this  report  in  terms  of  accuracy  and  utility.  To  the  degree  practical,  the  numerical  results  of 
the  evaluation  method  are  compared  against  theoretical  and  experimental  results.  The  intent  is  to 
demonstrate  the  achieved  performance  for  alternate  grid  resolutions  under  conditions  of  varying 
complexity. 


4.1  COMPARISON  AGAINST  A  BASELINE  Pss 

The  first  level  of  validation  under  consideration  is  a  comparison  of  calculated  values  for  the 
probability  of  successful  search  Pss  against  theoretical  values.  To  do  this,  a  “mow  the  lawn” 
type  of  search  pattern  is  used  where  the  search  trajectory  is  set  to  lie  entirely  within  the  search 
region.  In  this  example  the  search  speed  of  the  platform  is  set  constant  and  the  turn  rate 
C  =  S/Rd  is  constrained  to  achieve  constant  area  coverage  per  unit  time,  as  discussed  in 
appendix  C.  For  the  condition  of  uniformly  random  mine  placements  over  the  search  space  S, 
this  affords  a  direct  comparison  of  the  search  path  results  against  an  equivalent  long,  single-leg 
search  path.  The  predicted  Pss  yields  (as  a  function  of  time)  the  probability  of  finding  the  next 
occurrence  of  a  mine  within  the  executed  portion  of  the  search  path.  It  has  the  functional  form  of 
the  Koopman’s  random  search  law  (reference  3) 

Pss(t)  =  l-exp[-NMPD<|>(t)],  (107) 

where  NmPd  represents  the  expected  number  of  detected  mines  in  the  search  space  and 

(j,(t)  =  ^D+2RDS(t-tB)  ( 1 08) 

defines  the  area  swept  out  by  the  searcher  over  time  t  relative  to  the  size  of  the  search  space.  As 
the  actual  number  of  mines  is  unknown,  NmPdMO  yields  (as  a  function  of  time)  the  expected 
number  of  mines  detected  in  the  executed  portion  of  the  search  path. 

Figure  15  shows  the  sample  search  region  used  for  the  validation.  A  search  speed  of  10 
knots  is  used  with  a  search  time  of  8  hours  and  1 4  minutes.  Other  scenario  parameters  include 
an  update  rate  of  9  seconds,  a  sensor  detection  radius  Rp  =  300,  and  Pd  =  0.95,  yielding  an 
advance  ratio  of  Pa  =  6.  Using  a  value  of  Nm  =  2,  a  square  search  region  of  10,000  yards  per 
side  is  selected  to  produce  a  search  plan  that  is  not  dominated  by  straight  search  legs.  For  this 
experiment,  a  fixed  resolution  is  used  in  lieu  of  the  aforementioned  scaling  parameterization. 

Four  separate  numerical  evaluations  were  run,  using  spatial  grid  resolutions  of  Ax  =  { 100,  50, 

25,  12.5},  with  Ay  =  Ax  for  each  evaluation.  The  consequence  on  grid  cell  density  within  the 
sensor  detection  region  and  within  the  evaluation  region  is  illustrated  in  figure  16.  The  total 
number  of  grid  points  required,  that  is  Ntot  =  (2Nmax)2,  is  affected  accordingly  for  each  case. 
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Figure  15.  Search  Path  for  Baseline  Pss  Experiment 
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Figure  16.  Example  of  Grid  Spacing  Variation 


The  resulting  values  for  probability  of  successful  search  Pss  from  the  computational 
experiment  are  shown  in  figure  17.  All  four  grid  resolutions  are  depicted,  although  similarity  of 
results  makes  it  difficult  to  discern  them  in  the  figure.  The  red  dots  denote  the  calculated  values 
of  Pss  and  are  presented  at  a  sampling  of  every  50  time  updates.  They  are  superimposed  over  the 
blue  line  depicting  the  theoretical  value  of  Pss  using  equation  (107).  The  curve  realizes 
Koopman’s  random  search  law  as  described  in  section  2.3  for  the  case  where  the  expected 
number  of  mines  in  the  search  space,  E{Nm}  =  2,  is  known.  The  black  and  blue  dots  correspond 
to  the  case  where  the  actual  number  of  mines,  Nm  =  2,  is  known.  The  black  correspond  to  the 
probability  of  finding  both  mines,  while  the  blue  correspond  to  finding  either  of  them. 


Figure  1 7.  Probability  Calculation  Results 


Figure  18  shows  the  difference  error  between  the  Pss  calculation  and  the  theoretical  values. 
All  the  error  curves  show  the  effect  of  the  non-linear  mapping  from  area  coverage  to  probability. 
The  impact  of  the  platform  turns  in  executing  the  search  plan  is  least  discernible  (i.e.,  has  lowest 
error)  at  the  highest  resolution  (case  A)  and  becomes  more  pronounced  as  resolution  degrades. 

At  the  lower  resolutions,  the  impact  of  grid  cell  alignment  relative  to  the  search  plan  is  evident  as 
another  source  of  error.  In  the  lowest  resolution  (case  D)  where  the  update  evaluation  region  is 
smaller  than  the  grid  point  separation,  there  is  “jumpiness”  in  the  prediction.  Such  sizing 
conditions  do  not  present  a  theoretical  shortfall,  however,  as  the  probability  calculations  are 
essentially  set  functions  and  it  is  the  accuracy  of  the  aggregated  evaluation  regions  that  drive  the 
prediction  accuracy.  Each  of  these  phenomena  should  be  anticipated  in  selecting  grid  size  and 
scale  during  a  search  path  evaluation. 
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Figure  18.  Pss  Error  Over  the  Search  Path 


The  numerical  accuracy  of  the  four  cases  is  summarized  in  table  1 .  Here,  the  “Maximum 
Error”  refers  to  the  maximum  difference  error  evaluated  for  Pss  over  the  search  path  for  the 
given  grid  resolution.  The  “Coverage  Bound”  is  an  approximate  bound  reflecting  a  worst-case 
condition  that  is  provided  for  comparison.  This  bound  is  based  on  the  unusual  properties  of  this 
test  case.  Specifically,  with  turn  rate  C  =  S/Rn  and  an  aggregated  turn  angle  of  zero,  the  search 
trajectory  path  length  is  identical  to  that  of  the  two  coverage  bounds  on  either  side  of  the 
trajectory  and  that  of  an  equivalent  straight-line  path  had  no  turns  been  made.  To  show  this, 
observe  that  the  trajectory  path  length  over  the  turn  segment  is 

Lt  =6-  Rd  =  CATc  •  Rd  =SATc.  (109) 

As  the  outer  coverage  boundary  realizes  an  arc  of  radius  2Rd,  the  outer  boundary  path  length 
during  the  turn  is  twice  that  of  the  search  platform  trajectory,  while  the  inner  path  length  is  zero. 
Hence,  over  a  sequence  of  turn  and  countertum,  all  path  lengths  become  equal.  As  the  grid 
realizes  cell  boundaries  that  consist  of  parallel  lines  in  the  x  and  y  direction,  for  any  linear 
segment  of  the  search  path  the  maximum  number  of  grid  cell  interactions  for  a  search  boundary 
line  is  limited  by  N0  =  S  AT/R0.  For  curvilinear  motion  where  the  grid  resolution  is  much 

smaller  than  the  radius  of  curvature  of  the  turn,  the  crossing  of  the  boundary  path  across 
individual  grid  cells  is  approximately  linear  and  Ng  applies  as  an  approximate  limiting  condition. 
Additionally,  for  each  grid  cell  interaction  at  the  boundary,  the  maximum  error  in  the  coverage 
calculation  due  to  lateral  grid  overlap  is  on  the  order  of  one-half  of  the  cell  density  value,  with 
equivalent  sized  errors  distributed  equally  on  either  side  of  the  path  boundary.  Hence,  the 
approximate  bound  |<t>e(t)|  <  N(iAx:/a0  becomes  a  limiting  condition  on  the  aggregated  set  size 
error.  The  mapping  to  the  Pss  error  as  used  in  the  Coverage  Bound  is  thus 


44 
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|psse(t)|  =  |1  -  exp(-  NMPD<Kt))-  [l  -  exp(-  NMPD(<j>(t)  +  4>e(t»)]  | , 

=  |exp(-  NMPD<|>(t))[exp(-  NMPD<t>e(t))- 1]  | , 

<  exp(NMPH|<J)e(t)|)-l. 

In  table  1 ,  the  error  bound  is  significantly  larger  than  each  of  the  observed  maximum  errors. 
The  bound  is  meant  to  serve  as  a  worst-case  indicator  representative  of  the  situation  where  the 
search  trajectory  lines  up  in  the  worst  possible  arrangement  relative  to  the  underlying  grid.  More 
typically,  averaging  effects  lead  to  calculation  errors  accumulated  over  the  search  trajectory  that 
are  much  smaller. 


Table  1.  Pss  Baseline  Comparison  Results 


Case  A 

Case  B 

Case  C 

Case  D 

Multiplier 

1 

2 

3 

4 

Ax  /  Grid  Size 

12.5/640000 

25/160000 

50  / 40000 

100/10000 

Maximum  Error 

0.000043 

0.000193 

0.000712 

0.001638 

Coverage  Bound 

0.017703 

0.035406 

0.070812 

0.141623 

4.2  VALIDATION  VIA  SIMULATION  EXPERIMENTS 

This  section  presents  a  limited  set  of  simulation  results  to  authenticate  the  probability 
predictions  established  for  the  various  search  events  described  in  the  preceding  sections. 
Validation  is  achieved  over  a  Monte  Carlo  ensemble  of  repeated  trials  when  the  observed 
frequency  of  occurrence  approaches  the  predicted  probability  as  the  number  of  samples  in  the 
experiment  becomes  large.  These  experiments  are  intended  to  be  informative  as  well,  and  they 
serve  to  demonstrate  how  search  event  probabilities  must  be  interpreted  in  applying  the  various 
criteria. 

For  all  the  simulation  experiments,  a  grid  specification  of  {Nmax  =  400,  Nadv  =  3,  Pa  =  4}  is 
applied,  yielding  a  search  box  of  20000  yards  per  side  and  a  grid  of  640  thousand  cells,  with  a 
resolution  of  25  yards.  The  first  experiment  involves  a  simple  “mow  the  lawn”  search  plan  as 
above,  except  that  the  search  platform  broaches  the  search  region  boundary  before  turning  to  the 
next  leg.  This  enables  a  comprehensive  search  of  the  region,  as  illustrated  in  figure  19.  A 
uniform  mine  density  function  with  the  expected  number  of  mines  in  the  search  region  of  Nm  =  4 
is  applied,  with  a  constant  detection  probability  of  Pd  =  0.8  and  a  detection  radius  of  Rd  =  300 
yards.  The  numerically  computed  probability  assessments  are  shown  in  figure  20  for  a  search 
path  extending  to  three  sweeps  of  the  area. 


45 


Event  Probability 


Figure  19.  Search  Path  for  the  Uniform  Drop  Experiment 


Figure  20.  Search  Probabilities  for  the  Uniform  Drop  Experiment 
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First,  observe  in  figure  20  that  as  coverage  is  applied  outside  the  designated  search  area 
where  placement  likelihood  is  zero,  search  effectiveness  is  delayed  and  the  calculated  Pss  lies 
below  the  theoretical  value  conditioned  on  continuous  coverage.  Observe  also  that  the  next 
detection  probability  for  the  case  of  a  known  number  of  mines  in  the  search  region  lies  above  the 
Pss  curve.  This  precise  knowledge  of  the  tactical  situation  reduces  uncertainty  and  that  is 
reflected  in  the  prediction.  The  remaining  three  curves  relate  to  area  clearance  probability. 
Observe  that  the  spatial  average  clearance  probability  lies  consistently  above  the  probability  of 
the  known  number  of  mines  being  detected.  Similarly,  the  minimum  cell  clearance  probability  is 
also  shown. 

The  first  Monte  Carlo  experiment  focuses  on  next  detection  probability  and  its  dependency 
on  prior  knowledge  of  the  number  of  mines  Nm.  A  uniform  distribution  over  S  is  applied.  For 
the  known  number  of  mines  condition,  precisely  four  mines  are  placed  in  each  experimental 
realization.  For  the  unknown  number  case,  an  expected  value  of  four  mines  is  used  and  a 
Poisson  distributed  number  is  drawn  to  determine  the  actual  number  of  mines  to  place  at  random 
for  each  realization.  Mine  locations  are  tagged  to  grid  cells  and,  if  a  mine  is  encountered  during 
the  sequential  search  evaluation,  a  uniform  random  number  is  drawn  and  tested  against  a 
threshold  of  detection  probability  PD  to  determine  if  detection  occurs.  The  number  of  mines 
detected  over  the  search  path  is  recorded  and  a  test  statistic  is  set  to  1  if  a  detection  occurred  and 
to  0  otherwise.  A  partial  search  path  is  examined  for  this  experiment  with  results  generated 
when  the  predicted  Pss  attains  a  value  of  0.8.  This  yields  a  first  detection  probability  of  0.8725 
for  the  case  of  known  Nm  and  provides  a  context  where  the  predictions  are  distinct.  Figure  21 
presents  the  ensemble  average  of  this  test  statistic  as  a  function  of  the  number  of  ensemble 
realizations  (note  that  the  horizontal  axis  of  figure  2 1  is  logarithmic).  Observe  that  as  the 
number  of  realizations  approaches  1 000,  the  observed  frequencies  of  occurrence  converge  to 
their  predicted  values,  such  that  a  separation  between  outcomes  is  readily  apparent. 


Figure  21.  Partial  Search  Path  Mine  Detection  Results 


Next,  consider  a  clearance  search  path  comprising  an  entire  sweep  of  the  search  region.  For 
this  Monte  Carlo  experiment,  two  clearance  events  are  considered.  Given  precise  knowledge  of 
the  number  of  mines  Nm,  the  predicted  probability  of  detecting  all  of  them  is  compared  to  the 
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observed  frequency  of  occurrence.  To  do  this,  the  number  of  mines  that  are  detected  over  the 
search  path  is  recorded.  The  test  statistic  is  set  to  1  if  this  number  equals  the  known  value  and  to 
0  otherwise.  Results  are  plotted  in  figure  22.  The  other  clearance  event  involves  the  search 
effort  assessment  given  no  knowledge  of  mine  placement  within  the  search  region.  Specifically, 
the  frequency  of  occurrence  of  the  event  that  a  mine  exists  in  the  cell  and  is  not  detected  is 
recorded  for  each  grid  cell.  Hence,  both  placement  and  detection  events  require  recording,  with 
the  cell-specific  event  occurrence  given  by  the  ratio  of  placement  to  detection  events.  For  this 
experiment,  it  is  necessary  that  the  test  statistic  not  be  dependent  on  the  actual  number  of 
occurrences.  Otherwise,  it  would  depend  on  the  number  of  realizations  and  not  be  a  valid 
context  for  the  experiment.  The  cell-specific  occurrences  are  averaged  over  the  grid  and  the 
results  plotted  in  figure  22. 


Number  of  Samples 


Figure  22.  Clearance  Event  Frequency  of  Occurrence 


The  convergence  of  the  observed  experimental  frequencies  of  occurrence  to  their  predicted 
probabilities  is  apparent  in  figure  22.  The  achieved  tolerances  exhibited  in  figures  21  and  22 
remain  to  be  articulated,  however.  What  is  important  to  observe  is  the  difference  in  the 
performance  assessments  themselves.  The  search  effort  evaluation  presented  in  the  cell  average 
probability  corresponds  to  a  blind  search  where  it  is  unclear  what  the  actual  context  is  and, 
hence,  precisely  what  it  is  that  you  are  looking  for.  As  this  measure  is  developed  over  the 
probability  density  function  extracted  from  the  mine  density  ^-decomposition,  a  single  mine 
context  applies.  Thus,  the  result  provides  a  probability  that  if  a  single  mine  were  present  in  the 
region,  it  would  be  detected.  For  this  case  where  detection  probability  Pd  and  mine  density  X  are 
constant  over  the  region,  it  becomes  equivalent  to  the  sensor  characterization  independent  of 
geography.  Note  that  this  will  not  generally  be  the  case.  The  probability  of  detecting  all  mines 
in  a  search  region  presents  a  more  stringent  criterion  to  meet — it  is  contingent  on  having  that 
precise  knowledge  in  forming  the  evaluation  metric.  For  this  case,  the  search  is  presumably 
highly  directed  and  search  effort  applied  until  all  mines  have  been  detected  to  within  a  specified 
probability.  For  the  single  area  sweep  presented  in  figure  22,  the  predicted  probability  of  finding 
four  mines  is  0.4096.  This  contrasts  with  a  search  effort  threshold  of  0.8  that  was  applied.  In 
comparison,  a  Pss  next  detection  probability  calculation  using  an  expected  number  of  mines  of 
Nm  =  1  yields  a  predicted  value  of  0.5507. 
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The  remaining  experiments  are  conducted  using  non-uniform  mine  and  clutter  densities.  In 
addition,  a  field  structure  is  applied  in  developing  the  mine  density  function.  For  these 
experiments,  the  mine  field  consists  of  a  single  line  with  12  possible  mine  locations  along  the 
line.  The  expected  number  of  mines  populating  the  field  is  set  to  four.  A  random  mine 
placement  is  also  applied  with  an  expected  number  of  two,  yielding  a  total  expected  number  of 
six  mines  placed  in  the  search  region.  These  numbers  are  kept  artificially  low  so  as  to  extend  the 
partial  search  path  before  a  Pss  threshold  criterion  of  0.8  is  met,  as  the  intent  is  to  demonstrate 
the  properties  of  the  search  evaluation  capability.  The  partial  search  path  where  this  condition  is 
met  is  depicted  in  figure  23.  Also  shown  are  contours  corresponding  to  the  respective  mine  and 
clutter  density  functions. 

To  emulate  mine  classification  degradation  in  the  presence  of  clutter,  an  artificial  reduction 
in  detection  probability  in  the  form  of 

PD(x)  =  0.8-0.2-Mx)ACmax  (HO 

is  introduced  to  yield  an  unweighted  average  Pd  of  0.764  over  the  space  and  a  density  weighted 
average  of  0.732.  The  combination  of  spatial  variability  and  mine  field  structure  results  in  the 
predicted  probabilities  shown  in  figure  24.  The  impact  of  spatial  variability  on  the  next  detection 
probabilities  relative  to  that  of  the  baseline  Pss  (recall  that  the  baseline  is  conditioned  on  a 
uniform  distribution)  is  readily  apparent  for  this  partial  search  path.  The  separation  between  the 
next  detection  probabilities  assuming  a  known  or  unknown  number  of  mines  is  less  apparent  than 
in  figure  20.  This  is  due  to  the  increase  in  the  expected  number  of  total  mines  in  the  search 
space.  In  general,  the  next  detection  probability  curves  for  a  known  and  unknown  number  of 
mines  converge  as  the  expected  number  of  mines  Nm  increases. 


Figure  23.  Structured  Drop  Experiment  Mine  and  Clutter  Densities 
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Figure  24.  Structured  Drop  Experiment  Predicted  Probabilities 


Figure  25  illustrates  the  results  of  a  simulation  conducted  to  demonstrate  anticipated  mine 
detections  versus  detections  due  to  false  alarm  clutter.  The  clutter  density  spatial  Poisson 
process  is  colored  by  a  value  of  Pc  =  0.85,  denoting  the  probability  of  classifying  the  clutter 
properly.  The  unclassified  clutter  forms  an  independent  spatial  Poisson  process,  and  the 
detection  of  these  mine-like  objects  constitutes  a  false  alarm.  Figure  25  plots  the  ensemble 
average  of  the  number  of  detections  occurring  along  a  partial  search  path  that  is  terminated  when 
the  false  alarm  probability  reaches  a  value  of  0.8. 


Figure  25.  Detection  and  False  Alarm  Frequency  of  Occurrence 
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The  expected  values  of  the  number  of  mine  detections  and  false  alarms  encountered  along 
the  search  path  provide  additional  criteria  in  comparing  search  paths  under  geometrically  varying 
conditions.  Figure  26  shows  the  evolution  of  the  expectations  as  the  search  path  is  executed. 

The  end-case  ensemble  of  2000  simulation  runs  was  used  in  constructing  this  figure.  For  this 
depiction,  one  complete  sweep  of  the  search  region  is  executed.  The  figure  shows  the  calculated 
expected  values,  the  average  of  the  number  of  observed  detections  that  occur  for  each  time  step, 
and  a  set  of  bounds  corresponding  to  the  standard  deviation  of  the  ensemble.  For  this  example, 
the  results  demonstrate  the  dependency  on  spatial  variability  and  the  richness  of  the  false  alarm 
density  relative  to  that  of  the  mine  density. 


Figure  26.  Expected  Number  of  Detections  and  False  Alarms 
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5.  SUMMARY  AND  FUTURE  WORK 


This  report  has  detailed  the  development  of  a  search  path  evaluation  capability  for  objects 
that  may  be  placed  with  spatial  structure,  such  as  occurs  in  mine-hunting  applications.  The 
method  is  based  on  a  rigorous  extension  of  previous  work  on  distributed  search  theory,  whereby 
searcher  motion  is  combined  with  spatial  patterns  of  the  objects  of  search.  This  probabilistic 
view  of  the  spatio-temporal  search  process  is  consistent  with  uncertainty  representations  that  are 
characteristic  of  information  provided  for  the  assessment  of  realistic  mine-hunting  operations. 

The  method  is  appropriate  for  objects  occurring  randomly  or  in  spatial  patterns,  for  a  known  or 
unknown  number  of  objects,  and  with  or  without  clutter  contacts.  Furthermore,  the  method 
works  with  search  paths  of  arbitrary  complexity,  including  search  paths  that  self-intersect. 

In  the  process  of  developing  this  capability,  a  significant  degree  of  mathematical  rigor  was 
applied  to  create  the  event  spaces  over  which  probabilistic  calculations  are  based.  As  a 
consequence,  several  maps  were  developed  either  on  the  global  search  region  scale  or  on  the 
local  sensor  referenced  subspace.  These  include  the  respective  mine  and  clutter  density  maps, 
the  detection  and  non-detection  probability  maps,  and  the  respective  mine  field  maps  associated 
with  a  reference  placement  of  an  identified  mine.  These  maps  are  a  byproduct  of  the  search 
evaluation  procedure;  yet,  they  may  prove  to  be  useful  as  additional  information  to  an  operator. 
The  cases  of  prior  knowledge  on  the  number  of  mines  in  the  search  environment  or  its  expected 
value  were  considered.  When  mine  field  structure  is  assumed,  the  actual  or  expected  value  of  the 
number  of  mine  fields  is  applied  as  well.  A  variety  of  performance  results  were  developed  and 
presented  to  clarify  the  capability  in  terms  of  its  functionality,  efficacy,  and  accuracy. 

The  numerical  procedure  for  search  evaluation  employs  a  grid-based  realization  so  as  to 
readily  incorporate  environmental  influences  on  mine  placement  and  the  spatial  characterization 
of  clutter.  Additionally,  the  capability  integrates  a  degree  of  prior  knowledge  on  mine  field 
structure  into  the  evaluation  methodology.  This  modeled  structure  is  incorporated  into 
calculation  of  the  probability  of  finding  the  next  mine  for  the  cases  of  no  identified  mines  or  a 
spatially  diverse  set  of  known  identifications.  The  method  is  based  on  a  notion  of  search  effort 
swept  over  time,  so  it  is  directly  applicable  to  any  sensing  modality  that  can  be  represented  as  a 
probability  over  range  in  directions  relative  to  searcher  motion.  This  search  evaluation  capability 
does  not  address  the  fusion  of  prior  knowledge  of  mine  field  structure  with  the  location  of 
multiple  identified  mines  within  mine  field  proximity  to  each  other.  It  can  incorporate  such 
output  from  algorithms  that  do  perform  fusion  on  multiple  observations  if  it  is  provided  in  the 
prior  knowledge  specification. 

While  the  method  is  appropriate  for  robust  computation  of  problems  of  current  interest, 
there  are  a  number  of  extensions  that  can  improve  the  capability.  The  tessellation  of  the  search 
space  (to  represent  underlying  mine  and  clutter  density  functions)  was  limited  to  square  grids  for 
ease  of  development.  There  are  many  alternate  tessellations,  both  uniform  (such  as  hexagons) 
and  non-uniform  (such  as  different  sized  rectangles),  that  may  prove  to  be  computationally 
advantageous.  As  the  problem  scale  increases,  this  is  an  area  that  may  merit  consideration. 
Another  computational  complexity  issue  is  to  determine  the  appropriate  tradeoffs  between 
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effectiveness  of  computation  and  speed  of  computation,  which  can  be  computed  as  a  function  of 
grid  evaluation  parameters  using  available  Pareto  optimization  techniques  (reference  13). 

With  respect  to  the  search  path  input,  the  method  currently  assumes  a  perfectly  known 
search  path,  with  uniformly  sampled  time  intervals.  Obviously,  time  sampling  can  be  applied  to 
any  given  path  to  meet  the  uniform  sampling  intervals,  but  that  may  also  prove  to  be 
computationally  inefficient,  especially  for  search  paths  that  have  extended  periods  of  non¬ 
maneuvering.  Finally,  the  search  paths  are  assumed  to  be  known  precisely;  assessment  of  search 
effectiveness  under  the  anticipated  uncertainty  in  navigational  accuracy  of  the  search  platform  is 
of  critical  importance.  Since  the  search  effectiveness  is  based  on  probability  representations  of 
searcher  performance,  the  effects  of  positional  uncertainty  of  the  searcher  are  regularly 
incorporated  into  the  method;  this  is  a  subject  of  current  investigation. 

Future  extensions  of  this  search  evaluation  capability  include  the  development  of  additional 
probability  metrics  for  search  evaluation  beyond  probability  of  successful  search.  Of  primary 
importance  in  the  mine-hunting  application  is  the  representation  of  the  probability  of  obtaining 
queue  roots  that  provide  a  mine-free  (within  some  pre-determined  certainty)  path  of  given  width 
but  arbitrary  shape  between  two  points.  A  direct  extension  of  this  method  is  to  multiple 
searchers,  as  the  maps  that  underlie  the  search  computations  can  be  updated  simultaneously  for 
any  number  of  collaborating  searchers.  The  impact  of  data  fusion  between  such  searchers  is  also 
a  simple  extension,  given  a  known  description  of  the  impact  of  the  fusion  process  on  the 
corresponding  Pd  values. 

The  initial  uses  of  this  search  evaluation  capability  include  the  determination  of 
effectiveness  inside  of  various  path  optimization  strategies.  Such  applications  will  take 
advantage  of  the  robust  computational  nature  of  the  approach.  The  evaluation  of  various 
competing  fusion  strategies  for  multiple  searcher  mine  hunting  is  another  initial  area  of  use. 
Other  intended  uses  include  the  development  of  new  displays  that  represent  regularly  updated 
mine  location  likelihood  maps  as  search  progresses.  Within  the  operational  community,  the 
search  evaluation  capability  can  be  used  to  determine  the  impact  on  search  evaluation  obtained 
by  providing  more  (or  less)  information  about  the  expected  mine  field  characteristics.  All  of 
these  applications  are  under  consideration. 
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APPENDIX  A 

PROBABILITY  BACKGROUND 


A.l  REVIEW  OF  PROBABILITY  SPACES 

A  probability  space  represented  by  the  triple  (Q,  9,  P)  comprises  a  totally  finite  measure 
space ,  with  P  denoting  a  probability  measure  (references  8  and  14).  Here,  Q  denotes  the  set  of 
possibilities  under  consideration,  9  denotes  an  event  space  consisting  of  a  collection  of  selected 
subsets  of  Q,  and  the  probability  measure  P:  9  — >  [0,  1  ]  realizes  a  set  function  mapping  events 
E  e  9  to  numbers  denoting  the  likelihood  of  event  occurrence.  The  set  Q  is  often  described  as 
the  set  of  possible  outcomes  of  some  experiment  where  Q  is  exclusive  and  exhaustive;  that  is, 
outcomes  can  consist  of  only  one  member  of  the  set  and  no  other  outcomes  of  the  experiment  are 
possible.  In  contrast,  events  often  consist  of  a  set  of  outcomes  and  the  event  space  may,  or  may 
not,  include  individual  outcomes  as  singleton  events. 

For  9  to  be  a  valid  event  space,  9  must  comprise  a  a-algebra  and  every  element  of  9  must 
be  P-measurable.  For  any  collection  9  cz  5*(Q),  where  5*(Q)  represents  the  power  set  of  all 
possible  subsets  of  Q,  9  is  a  a-algebra  if  and  only  if: 


(i)  <J>  e  9. 


(ii)  If  E  e  9,  then  ~E  e  9. 


(iii)  For  any  countable  sequence  (Ej^  of  elements  of  9, 


ac 


i=l 


In  the  above,  ~E  denotes  the  set  complement  of  E  with  respect  to  Q,  and  <|>  is  the  empty  set.  As  <J> 
must  be  in  the  event  space,  so  must  be  its  complement  Q.  Also,  from  (ii)  and  (iii),  all  countable 
intersections  of  events  must  be  in  the  event  space.  If  9  is  a  valid  a-algebra,  then  (Q,  9)  is  a 
measurable  space.  For  P  to  be  a  valid  measure  on  (Q,  9),  the  following  must  hold: 


(i)  P(< \>)  =  0. 


(ii)  />(E)  >  0  for  all  E  e  9  . 


(iii)  For  any  countable  sequence  (E,)jH  of  elements  of  9  where  the  Ej  are  disjoint. 


If  these  conditions  are  met,  the  triple  (Q,  9 ,  P)  is  called  a  measure  space.  The  measure  space 
becomes  a  probability  space  when  P  (Cl)  =  1 ,  with  P  being  a  probability  measure. 


A-l 


Observe  that  the  smallest  possible  event  space  (i.e.,  a-algebra)  consists  of  the  set  {<)),  Q}. 
However,  this  a-algebra  is  non-informative  and  no  probability  questions  regarding  a  particular 
subset  of  Q  can  be  asked.  What  is  important  for  purposes  of  search  path  evaluation  is  that  event 
spaces  are  constructed  that  are  specific  to  the  set  of  probability  questions  to  be  answered.  That 
is,  given  a  probability  question  (e.g.,  the  probability  of  finding  a  mine  along  a  given  search  path 
within  a  specific  interval  of  time),  an  event  space  &  is  constructed  that  yields  the  smallest 
a-algebra  necessary  to  answer  the  question. 


A. 2  REVIEW  OF  SPATIAL  POISSON  PROCESSES 


The  spatial  Poisson  process  is  a  multi-dimensional  extension  of  the  set  of  arrival  times 
associated  with  a  one-dimensional  Poisson  process.  For  our  purposes,  the  arrival  characteristic 
concerns  the  placement  of  objects  in  a  two-dimensional  space.  Objects  persist  in  the  space,  so 
the  experimental  outcome  occurs  only  once.  However,  given  any  arbitrary  region  AcS,  the 
number  of  objects  occurring  in  that  region  of  the  search  space  N(A)  becomes  a  Poisson 
distributed  random  variable.  The  details  of  this  process  are  given  below. 

Let  A  be  a  collection  of  Lebesgue  measurable*  regions  {A} a  e  a  within  the  search  space  S. 
Then,  {N(A)}a  £  a  is  a  homogeneous  spatial  Poisson  point  process  with  constant  intensity  X  if: 


(i)  For  each  A  e  A,  N(A)  ~  Poisson(X|A|). 

(ii)  For  every  finite  collection  {A|,  A2,...,  An}  of  disjoint  subsets  of  S, 

N(Aj),  N(A2),...,  N(An)  are  independent  random  variables. 

There  are  alternate  axiomatic  conditions  that  the  spatial  Poisson  process  satisfies.  These  are: 


(0 


If  A|,  A2,  ...,  An  are  disjoint  regions,  then  N(A|),  N(A2),  ...,  N(An)  are  independent 
random  variables  and 

f  n  \  n 


N 


Ik  =Zn(a,). 


i=l  7  i=l 


(ii)  The  probability  distribution  of  N(A)  depends  on  set  A  only  through  its  size  |A|. 

(iii)  There  is  zero  probability  of  points  overlapping: 

Pr(N(A)  >  1)  x 
il™  Pr(N(A)  =  1) 


Then,  for  any  spatial  region  AcS,  the  probability  of  finding  k  objects  within  A  for 
k  =  0,  1,  2,  ...,  is 


*Any  countable  union  of  intervals  in  P,  or  rectangles  in  P:,  is  Lebesgue  measurable. 
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Pr(N(A)  =  k)  =  ^-exp(-X|A|). 

Next,  the  requirement  of  using  a  Lebesgue  measure  on  P2  is  relaxed  and  variability  of  the 
intensity  parameter  X,  over  the  search  space  S  is  allowed.  Let 

A(A)  =  J^X(x)dx 

for  all  Lebesgue  measurable  AcS.  Denote  the  set  function  A(A)  as  the  mean  measure  of  A  and 
variable  X  as  a  characterization  of  a  non-homogeneous  spatial  Poisson  process.* 

The  axiomatic  requirements  for  such  a  process  become: 

(i)  For  each  A  e  A,  N(A)  ~  Poisson(A(A)). 

(ii)  For  every  finite  collection  {A|,  A2,...,  An}  of  disjoint  subsets  of  S, 

N(A|),  N(A2),...,  N(An)  are  independent  random  variables. 

Two  important  theorems  regarding  spatial  Poisson  processes  of  particular  relevance  to 
search  path  evaluation  are  stated  here  in  terse  form  and  without  proof.  Readers  can  refer  to 
reference  8  for  a  more  complete  coverage  of  the  topics. 

Superposition  Theorem:  Let  n  1  and  n2  be  independent  spatial  Poisson  processes  with 
intensity  functions  X|  and  X2.  Then,  the  set  n  =  n  1  n  n2  is  a  spatial  Poisson  process  with 
intensity  function  X  =  X|  +  X2. 

Coloring  Theorem:  Let  n  be  a  non-homogeneous  spatial  Poisson  process  with  intensity 
function  X(x).  Let  each  point  of  n  be  “colored”  independently  as  green  at  x  with  probability  y(x) 
and  red  with  probability  p(x)  =  1  -  y(x).  Let  T  and  P  be  the  set  of  points  that  are  colored  green 
and  red.  Then,  T  and  P  are  independent  spatial  Poisson  processes  with  respective  intensity 
functions  y(x)X(x)  and  p(x)X(x). 

These  properties  and  theorems  allow  the  construction  of  mine  density  functions  and  sensory 
characterizations  that  vary  over  the  search  space.  To  show  this,  let  X(x)  be  a  mine  density 
function  (i.e.,  intensity)  associated  with  a  spatial  Poisson  process  n.  Then,  for  each  point  of  n, 
let  Up  be  the  occurrences  colored  by  the  detection  event,  with  probability  Pd(x),  and  let  IIq  be 
the  occurrences  colored  by  non-detection,  with  probability  Qd(x)  =  1  -  Pd(x).  Then,  from  the 
above  theorems,  rip  and  TIq  constitute  independent  spatial  Poisson  processes  with  respective 
density  functions  PD(x)X(x)  and  ( 1  -  PD(x))X(x). 


*Note  that  when  X  is  constant,  A(A)  =  j\Xdx  =  X  j\dx  =  X|  A  | . 


A-3  (A-4  blank) 


APPENDIX  B 

JOINT  DENSITY  DEVELOPMENT  FOR  GAUSSIAN  FIELDS 


This  appendix  provides  the  development  for  the  joint  probability  density  function  discussed 
in  section  2.2  for  the  case  of  a  field  structure  comprising  Gaussian  components.  Each  structure 
component  is  modeled  as  a  normally  distributed  random  variable  Xj,  j  =  1, Nm  placed  about 
the  common  reference  point  xref  in  the  search  space  S  <z  R2.  Let  ^  =  xref  +  AXj  designate  the 

intended  position  of  object  mj  with  displacement  AXj  from  the  reference.  Let  each  reference 
conditional  placement  be  distributed  as  Xj  ~  N[Xj,I2a2]  and  assume  that  tail  probability  outside 
S  is  negligible.  As  shown  in  equation  (49),  the  joint  conditional  density  takes  the  form 


f(x„...,xNJ  xref)  =  nfJ(xJ|xref), 


j=l 


=  [2™p]NMexP  -T^r5((xj  +(yJ-yJ)2) 


2ct;  j=i 


(B-l) 


Let  the  field  placement  density  also  be  Gaussian,  xref  ~  N[x0,I2ar2ef  ] .  Then,  the  unconditional 

joint  density  can  be  developed  via  marginalization  as  in  equation  (41).  This  development  is 
tedious,  however,  and  it  is  more  insightful  to  map  the  conditional  density  into  a  form  where  the 
unconditional  joint  density  is  readily  recognized,  as  in 


f(x,,...,xNM,xref)  =  f(xref  |x,,...,xNM)-f(x,,...,xNM), 


(B-2) 


Here,  f(  xref  |  x,,...,xN  )  represents  the  posterior  probability  density  function  of  where  the 

reference  point  may  be  located  given  the  field  placement  prior  and  independent  observations  of 
the  Nm  distinct  objects  in  the  structure.  The  notation  -  Xj  -  Ax}  =  xref  +  n  t  is  used  to  denote 

observations  on  reference  location,  with  ~  N[0,I2a2]  as  above.  Let 

Z  =  [z,,z2,...,zNJ'  =  A-xref +  [n,,n2,...,nNJT  (B-3) 


denote  the  vector  of  such  observations  with  A  =  [  I2, 12,...,  I2]7  and  covariance  matrix 
Rz  =  a2  •  I2N^  .  Then,  the  reference  conditional  and  prior  densities  take  the  form 


f (x, , . . . ,  xN  I  x ref )  = 


tor 


expl-itZ-Ax^rR^Z-Ax,.,] 

V  ^ 


(B-4) 


f(xref)=,  2  exp 


2na 


ref 


f  1  *\ 

--[x0-xrefr  Ro'[x0-xref] 


(B-5) 
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where  R0  =  ajef  •  I2 .  The  joint  density,  including  field  and  object  placement,  becomes 


f(x,,...,x  ,xref)  = 


1 


[2Traj]NM2KCT;f 


—exp 


^  [z  Axre(  J  Rz  [Z  Axref]  ^  [x0  xref]  R0  [x0  xrd  ] 


Let  L  denote  the  exponent  of  equation  (B-6)  without  the  scalar  term,  -1/2.  Then, 

L  =  [Z  -  Axref  r  R-'[Z  -  Axref  ]+  [x0  -  xre  J  R0'  [x0  -  xref  ], 

=  x^[R0-1+ATR>]xref-2x^f[R-,x0  +  ATRz-,z]+ZTR;,Z  +  xjR-|x0. 

With  pp  =  ap/arcf  and  following  the  standard  development  in  reference  15,  the  notation 


P*/z  =  R<>  +  A  Rz  A, 


U*  1  .  _  NMa;f+a^  ,  NM+p;  , 

2  _ 2  l2  ~  7  1 


—  I2  +  I— l 12  = 


;  ref 


na: 


°ref°p 


a: 


"ref 


=  [Ro1  +  atr;'a]"'  [r-'x0  +  atr;'z], 


— x0  +  I-rzj 


NM+Pn 


'  ref 


h  a: 


■xn  +  - 


N, 


xi  2  0  vr 

nm+Pp  Nm+P 


pLNm  h 


S(x,-Ax,) 


is  used  to  denote  the  respective  information  matrix  and  conditional  mean  of  the  reference 
posterior  density  f( xref  |  x,,...,xN  ).  Completing  the  square  in  equation  (B-7)  yields 


L  =  x‘ fp;'z  xref  -  2xrefP“'zxrcf  +  z  R  'Z  +  XqRq  x0, 

=  [Xref  —  Xref  Px/Z  [x ref  —  Xref  ]+  Z  Rz  Z  +  X()  R0  XQ  —  Xref  Px/zXref . 

That  part  of  L  that  does  not  correspond  to  the  reference  point  posterior  becomes 


Lz  =  Z  R  Z  +  x0R-  x0  -  xret-P;/zxref, 


=  ZTR;'Z  +  XpR0'x0  -  [r-'x0  +  AtR”'z]'  Px/z[ro'x()  +  A'  Rj'z] 

=  zt[r-,-r>px/zatr-,]z-2ztr;1apx/zRo,Xo  +  x:[r-,-r-,px/zr0i]x0. 


To  complete  the  square  in  Z,  let 


P- ■' =  R-1  -  R-'APX  zatr:', 


k  x/zx 

±1  1  «  '  1.1  - 1 

O'  'Nm  NMCT^f  +cr  CT  Op 

p  p  M  rei  p  p  p 


'  ref 


2N 


M  NMa;f  +  a- 


-AAr 
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(B-6) 


(B-7) 


(B-8) 


(B-9) 


(B-10) 


(B-ll) 


(B-12) 


denote  the  information  matrix  of  the  marginalized  density.  Covariance  matrix  Pz  is  readily 
obtained  by  applying  the  matrix  inversion  lemma  to  equation  (B-I2)  and  simplifying  terms. 
Doing  so,  one  obtains 


P,  =  R,  +  a[p-',  -  atr;'a]'1  at  =  r,  +  ar„at, 

=  +  a>AT. 


Using  this  form,  the  marginalized  density  mean  value  simplifies  accordingly  to 

z  =  [r;1  -  r;'apx/zatr;,]"'[r;'apx/zr'ix0]. 


=  kl2NM+<4AAT] 


1  CTp  1  Ax 

CTp  ^  M  +  Pp  aref  °’ 


-A  +  - 


1 


NM+P  Nm  +p: 


•aata 


x0  =  Ax0, 


(B-14) 


where  the  substitution  ATA  =  NMI2  has  been  used.  Completing  the  square  yields 


L,  =  [z  -  if  p;1  [z  -  i}+  xJ[r-'  -  R0-'p,„R,']x„  -  zTp;'z. 

Rearranging  equation  (B-8)  and  substituting  A 1  Rz'  A  =  Pjz  -  R0' , 

x»[r»'  -R«'p,/iR»'K-ztp;'z, 

RJ1  -  atr;'a  -  Ro'p.^r;1  +  a  r  'ap,  a  r,'  a]x„, 
r-'  -  a'r,'a  -  r-'p.^r-1  +  [p;,1,  -  r;']p.,.atr;1  a]x#, 
R,"-R,"p..R;,-R.'p.„air;ia]x,„ 

Ri'-Ri'P.!kl  +  ATR;1  a]x„ 

R0  —  R0  ]x0, 


(B- 15) 


=  x 
=  x 
=  x 
=  x 

=  X0 

=  0. 


(B-16) 


Hence, 

Lz=[Z-Ax0J,Pz-,[Z-Ax0]. 


(B-17) 


To  return  the  exponent  to  an  expression  in  the  respective  Xj,  recall  that  z-  =  x ^  -  Axj .  Let 
the  vector  of  object  positions  and  mean  values  be  defined  as 


X  —  [ X| , Xj, . . . , X|^^  ]  ,  X  -  [x0  +  Ax,,x0  +  Ax2,...,x0  +  AxNm  ]  . 


(B-18) 
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Then,  the  quadratic  form  in  equation  (B-17)  is  equivalent  to 


l,  =  [x-xf  p;'[x-x]. 


(B- 19) 


It  remains  to  show  that  the  normalizing  factors  associated  with  the  reference  conditional  and 
marginalized  densities  of  equation  (B-2)  produce  that  of  equation  (B-6).  This  is  achieved  when 

2ic|Pj*  [27t]NM|Pz|T  =  [27tCTj]NM27tCT^f .  (B-20) 

From  equation  (B-8),  the  determinant  of  the  reference  posterior  covariance  matrix  becomes 


P 

rx/z 


NM  +P 


nm+p; 


( B-2 1 ) 


p  J 


The  marginal  covariance  matrix  Pz  of  equation  (B-13)  is  decoupled  in  the  x  and  y  Cartesian 
coordinates.  Hence,  if  Pz  is  rearranged  into  block  diagonal  components  corresponding  to  the  two 
dimensions,  then 


P/  0 

lp2l= 

zx 

=  lp  1 

p 

0  pZy 

1  z»l 

zy 

(B-22) 


Each  coordinate  covariance  matrix  takes  the  form  of  a  symmetric  compound  matrix  of  the  form 
pz,  =  Pz  =ap!NM  +wwT ,  (B-23) 


where  w  =  [aref,aref,...,aref]T  is  of  size  (NM  x  1).  Hence,  as  in  reference  15, 


P,  = 


-  cr 


2Nv 


f  XT  2  ^ 

)+  NM°ref 


V 


=  a:N"a2 


^Pp  +  NmA 


ref 


P  / 


(B-24) 


p  / 


Substituting  equations  (B-21),  (B-22),  and  (B-24)  into  equation  (B-20)  verifies  the  equality. 
Thus,  the  marginal  density  for  object  placement  takes  the  form 


f(x,,..„xNJ  = 


[2ti]Nm 


exp 


:-x]'p;' 


x] 


(B-25) 
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Rewriting  the  exponential  in  equation  (B-25)  with  the  form  of  Pz  1  from  equation  (B-12)  leads  to 
the  following  final  form  for  the  joint  density: 


f(x„...,xNM)  = 


[2™p]Nm 


\ 


Nm+P 


xexp 


-  M 
( 


expj- 


p  / 


2<jp 


I(xJ-xJ)2+(yJ-yj)2 

U=> 


2a; 


1 


+  Pp  J 


fNj,  „  V  fNj) 

Ix-x,  +  Iy,-y, 

J=>  J  V  J=> 


(B-26) 
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APPENDIX  C 

GOAL-ORIENTED  CURVILINEAR  MOTION 


This  appendix  describes  the  goal-oriented  curvilinear  motion  (GOCM)  model  that  is 
integrated  as  an  external  process  within  the  search  evaluation  capability.  The  essence  of  GOCM 
kinematic  modeling  is  as  follows.  All  GOCM  vehicular  motion  is  modeled  as  being  achieved 
under  the  condition  of  constant  angular  and  radial  acceleration.  That  is,  let  C  and  S  denote 
constant  change  of  course  and  speed  rates,  where  C  and  S  represent  the  vehicle  course  and  speed, 
respectively,  as  functions  of  time.  Then,  for  an  initial  time  to  and  arbitrary  time  t, 

C(t)  =  C(t0)  +  (t-t0)-C  and  S(t)  =  S(t0)  +  (t-t0)S.  Using  x(t)  =  [x(t)  y(t)  x(t)  y(t)f 

to  denote  a  Cartesian  coordinate  representation  of  platform  position  and  velocity,  then 

C(t)  =  tan  1  (x(t)/y(t))  and  S(t)  =  -y/x(t)2  +  y(t)2  .  To  achieve  constant  velocity,  course  and 

speed  rates  are  set  to  zero.  This  condition  is  indicative  of  a  goal  state  satisfaction.  More 
generally,  the  desired  course  and  speed,  CGoai  and  SGoai,  will  deviate  from  current  values  and 
induce  a  maneuvering  trajectory.  Consider  AT  as  a  search  evaluation  update  interval.  Let  AT i 
denote  the  time  necessary  to  achieve  the  first  goal,  regardless  of  which  goal  parameter  it 
corresponds  to.  With  maneuver  times  defined  as 


Min 

f|Ca.,-C(t„)|  J 

,  ATS  =  Min 

MSc,  -s(t„)|  ATj 

l  c  ’  J 

Is  ’  J 

and  with  AT,  =  Min(ATc,  ATS),  the  GOCM  equations  of  state  can  be  summarized  as 


(C-l) 


x(t)  =  <T>(t0,  t)  x(t0). 


1 

0 

A  AT 

BAT 

0 

1 

-BAT 

A -AT 

0 

0 

Dc 

Ds 

0 

0 

-Ds 

Dc 

(C-2) 


where  Dc  =  (1  +  yATs)cos(CATc) ,  Ds  =  (1  +  yATs)sin(CATc)  ,  y  =  S/S(t0) , 
8  =  0  +  yATs )  •  (AT  -  ATC )  - y  •  (ATS  -  AT,)2/2 ,  and 


a  at  /.  ,T,sin(CATc) 
AAT  =  (l  +  yAT,) — y  c/-y 


BAT  = 


1  -  (1  +  yAT, )  cos(CATc ) 


1  -  cos(C  AT, ) 
C2 

sin(CAT,) 


+  8cos(CATc), 
+  8sin(CATc). 


(C-3) 


The  GOCM  model  will  produce  search  trajectories  that  exhibit  a  high  degree  of  regularity. 
In  particular,  the  sensor  coverage  methodology  discussed  in  section  3  is  applicable  if  the  update 
interval  AT  is  selected  to  constrain  course  and  speed  changes  to  nominal  values  where  trajectory 
segments  are  approximately  linear.  Platform  turn  rates  can  be  set  to  achieve  specific  search 


C-l 


patterns  such  as  a  “mow  the  lawn”  path  plan.  Specifically,  at  the  completion  of  a  search  leg, 
maintaining  the  speed  S  and  setting  CgobI  to  C(t)  +  180°  and  C  =  S/RD  will  produce  adjacent 
swaths  of  width  2RD.  In  addition,  the  area  covered  by  the  search  swath  from  such  a  trajectory  is 
identical  to  that  obtained  had  a  straight-line  path  been  executed.  To  show  this,  observe  that  the 
area  covered  by  a  sensor  of  detection  radius  Rd  during  a  turn  of  0  radians  following  a  trajectory 
with  radius  of  curvature  of  Rc,  Rc  >  Rd  is 

At  =e/2-((Rc  +  RD)2-(Rc-RD)2)=e-2RcRD.  (C-4) 

A  straight-line  trajectory  of  swath  width  2Rp  will  cover  an  area  of 

As  =  2  R  D  •  SATC  .  (C-5) 

With  Rc  =  Rd  and  ATC  =  0/C  =6  R„/S , 

At  =20R2r  =2Rd  SATc.  (C-6) 

With  speed  S  constant,  the  constraint  C  =  S/Rr)  yields  a  constant  area  coverage  per  unit 
time  for  non-overlapping  search  trajectories  regardless  of  how  often  directed  turns  are  made.  As 
the  probability  assessments  of  section  2  are  measures  on  sets  relatable  to  area  coverage,  this 
result  is  useful  for  validating  the  calculated  search  performance  against  theoretical  predictions 
using  constant  velocity  search  paths.  Search  path  length  LT  =  S  •  ATC  =  0  RDis  also  identical  for 
the  straight  and  curvilinear  trajectories  under  these  conditions. 
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